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1. Introduction. The impedance concept is the foundation of engineering trans- 
mission theory. If wave guides are to be fully utilized as transmission systems or 
parts thereof, their properties must be expressed in terms of appropriately chosen 
impedances or else a new transmission theory must be developed. The gradual ex- 
tension of the concept has necessitated a broader point of view without which an 
exploitation of its full potentialities would be impossible. 

In the course of various private discussions, I have found that there exists some 
uneasiness with regard to the applicability of the concept at very high frequencies. In 
part this may be attributed to relative unfamiliarity with the wave guide phenomena 
and in part to the evolution of the concept itself. Some particular aspects of the con- 
cept have to be sacrificed in the process of generalization and although these aspects 
may be logically unimportant, they frequently become psychological obstacles to 
understanding in the early stages of the development. For this reason I am going to 
devote several sections of this paper to a general discussion of the impedance concept 
before passing to more specific applications; then by way of illustration I shall prove 
that an infinitely thin perfectly conducting iris between two different wave guides 
behaves as if between the admittances of its faces there existed an ideal transformer. 
This theorem is a generalization of another theorem which I proved several years 
ago to the effect that when the two wave guides are alike, the iris behaves as a shunt 
reactor. Actual calculation of the admittances and the transformer ratio depends on 
the solution of an appropriate boundary value problem. 

More generally, wave guide discontinuities are representable by 7-networks. In 
some special cases these networks lack series branches and in other cases, the shunt 
branch. 

2. Evolution of concepts. Concepts evolve. It is a long way from the primitive 
to the modern number concept. The primitive number was an integer, a concrete 
integer at that. In some primitive languages there is no word corresponding to “two.” 
There are words meaning “two men,” “two horses,” etc.; but the concept of “two” 
applying either to men, or to horses, is lacking. To a primitive mind the difference 
between a class comprised of two men and a class comprised of two horses over- 
shadowed the similarity. Seeing similarities requires a degree of abstraction. A re- 
sistance to abstract ideas seems to be a characteristic of human minds even in modern 
times; only the modern mind is quicker to overcome it. An example, pertinent to 
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our present discussion, is the following excerpt from a paper Derivation and discussion 
of the general solution of the current flowing in a circuit containing resistance, self- 
inductance and capacity with any impressed electromotive force, by Frederick Bedell and 
Albert C. Crehore, published in the Journal A.I.E.E., 9, 340 (1892): 

“From the analogy of this equation to Ohm’s law, we see that the expression 

[R?+(1/Cw—Lw)?]'/? is of the nature of a resistance, and is the apparent resist- 

ance of a circuit containing resistance, self-inductance and capacity. This expression 

would quite properly be called ‘impedance’ but the term impedance has for sev- 
eral years been used as a name for the expression [R?+L*w?]'/?, which is the 
apparent resistance of a circuit containing resistance and self-inductance only. 

We suggest, therefore, that the word ‘impediment’ be adopted as a name for the 

expression [R?+(1/Cw— Lw)?]"/? which is the apparent resistance of a circuit con- 

taining resistance, self-inductance and capacity, and the term impedance be re- 

tained in the more limited meaning it has come to have, that is, [R?+L%w?]'/2, 

the apparent resistance of a circuit containing resistance and self-inductance only.” 
The name “impediment” was not adopted. Apparently, it was soon understood that 
if one really wished to emphasize the difference between the impedances of various 
circuits, one could simply describe the circuits and, therefore, for most purposes, it 
was best to emphasize the similarity rather than the difference. And only ten years 
ago there were some who objected to the use of the word impedance for the ratio 
E/H in an electromagnetic wave and who wanted a new word for it. 

The word “number” now includes fractions, negative numbers, irrational numbers 
and complex numbers; the impedance is now a complex number, and not its absolute 
value as originally intended. There are mechanical impedaices, acoustic impedances, 
electromechanical impedances, and finally impedances associated with any wave no 
matter what its physical nature happens to be. The impedance is now the force/re- 
sponse ratio when the force and response are harmonic functions of time and are 
represented by complex exponentials. Around this concept has grown the transmis- 
sion theory of force and response in linear systems. The principal tool of this theory is 
the theory of functions of a complex variable. This theory is used for engineering 
purposes as in the design of filters, equalizers, and other transmission systems with 
prescribed desired properties; and with equal advantage it may be used for general 
transmission studies. In this paper I am particularly concerned with fundamental 
ideas applied to wave guides and wave guide elements. . 

3. General discussion of impedance and admittance. Superficially, it may seem 
that the impedance concept does not apply to wave guides or if it does it is quite 
different from the concept as applied to ordinary transmission lines. Actually there is 
no significant difference; whatever difference there exists is largely psychological 
rather than logical. In wave guides a characteristic impedance has to be associated 
with each transmission mode. At first the existence of various transmission modes 
may strike one as a feature which distinguishes high frequency wave guides from low 
frequency “ordinary” transmission lines; but soon one will realize that even in ordi- 
nary transmission lines it is usual to distinguish between different modes of trans- 
mission. Consider, for instance, parallel wires at the same height above ground; there 
are two obvious transmission modes recognized by communication engineers; in one 
the currents in the wires are equal and flow in opposite directions and in the other 
they are equal, flow in the same direction and return through ground. It is the exist- 
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ence of these two transmission modes that accounts for the important engineering 
difference between balanced and unbalanced transmission lines. Similarly, there are 
two obvious transmission modes in a shielded parallel pair. The field patterns are 
different for different transmission modes and the characteristic impedances are usu- 
ally different. The existence of transmission modes is not peculiar to hollow tubes 
and other structures which have become prominent in high frequency transmission; 
high frequency transmission studies make us merely aware of the fact that any 
physical wave guide, whether a coaxial pair or a shielded pair or a hollow tube, 
admits of an infinite number of transmission modes with their characteristic field 
patterns, characteristic impedances, and propagation constants. 

Another cause of worry to some is a degree of indeterminacy connected with an 
impedance and its associated quantities. The characteristic impedance of a wave 
guide may be defined in a number of ways giving different values. For each oscilla- 
tion mode a cavity resonator behaves as an ordinary circuit comprised of inductance 
and capacitance; but different values are obtained, depending on how L and C are 
defined. This indeterminacy is really inherent in these conceptions but in elementary 
theory it is not stressed for the simple reason that no occasion arises for such stressing. 
In the final analysis, this indeterminacy is of the same kind as that involved in the 
essential arbitrariness of units and is related to the fact that properties of analytic 
functions are not affected by a constant factor. Putting it in the language of transmis- 
sion theory, the essential properties of impedance functions are not affected by ideal 
transformers. If we have a closed box containing an electric network with two acces- 
sible terminals and if we measure a resistance R across these terminals, we cannot 
be certain that the box contains a resistance R; it may contain a resistance R/10 
which is then boosted to R by an ideal transformer. It does not really matter which 
is the case. Similarly, if the measurement seems to indicate that in the box we have 
a tuned circuit with an inductance L in series with a capacitance C, we may actually 
have a tuned circuit with an inductance 3Z and the capacitance 2C in the secondary 
of an ideal transformer which then doubles the impedance. More generally, the im- 
pedance function is defined by its zeros, infinities, and other singularities except for a 
constant. 

If V is the voltage across an impedance Z, J the current through Z, and W the 


complex power, then 
Vv* 


07* ’ (3-1) 





V = ZI, W = 3VI* = 32ZII*, W= 


where the asterisk is used to designate conjugate complex numbers. Now suppose 
that our voltmeters and ammeters contain concealed ideal transformers; then “Z” 
will have different values in the above equations and we shall have 


VvV* 
2Z% Vv 





V=Zyixl, We=will*, W= (3-2) 


These new equations are in effect various definitions of impedance and admittance 
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? V 2W ; VV* 
ihe a oe (3-3) 
; 1 ; 1 1 sin 
Tie Vyas =——>? Ywv = >—° 

Zyit Zw Zwv 


Ordinarily, we make sure that there are no concealed ideal transformers in our 
measuring instruments. Furthermore, at low frequencies we can measure the voltage 
across the total capacitance! and the total current through the inductance. There 
seems to be no question about the meaning of “V” and “J” and it so happens that in 
this case we are led to equations (3-1). However, in a section of a transmission line 
or in a cavity resonator the capacitance and inductance are not localized and we are 
forced to recognize the existence of a certain amount of indeterminacy. There is no 
harm in this indeterminacy; it does not really matter in which of the following two 
forms we decide to write the expression for power 

if 
W = 3ZI1I*, or W = 3(nZ)— —» (3-4) 
nn 
so long as we know how to compute it. 

Just as ideal transformers in our “ammeters” and “voltmeters” transform equa- 
tions (3-1) into equations (3-2) in the case of “ordinary” networks, they may be used 
to transform equations (3-2) into (3-1) in the case of wave guides and networks with 
distributed constants. 

4. General impedance relations. Eliminating V, J, and W from (3-2), we have 
the following equation connecting various impedances 


Seales _ Se Aes. (4-1) 

If the impedances are real, then 
Zwizwy = Zia (4-2) 
In equations (3-2) V and J may be arbitrarily chosen values of the voltage and 


current associated with a given impedor. If we choose a given definition for J, we can 
define a voltage 

Vw = Zyl (4-3) 
for which equations (3-1) will hold and the impedance Zyw,r will become the only 
impedance associated with the impedor. We can also define 


V 
’ (4-4) 





Iw.v = o 
“WV 
so that again we shall have equations (3-1) with Zw,y as the sole impedance. 
Since the power is an invariant we have 
_ * 
J wi Iw.v 


Vust = View Se (4-5) 
V r” 





1 Or almost across the total capacitance. 
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From (4-3), (4-4), and (4-5) we have 


| Vw | Iw.v | Zw. Yww 
J Ley Zw Ywu 
It is now evident that we can base our calculations on any particular voltage-current 


pair and then, whenever desirable, we may pass to any other pair simply by inserting 
in our transmission diagrams an ideal transformer with a proper impedance trans- 








formation ratio. 

5. Characteristic impedances and admittances of wave guides. The basic imped- 
ance associated with the mth transmission mode in a wave guide is defined as the ratio 
of the transverse electric to the transverse magnetic intensity 

K.« Fen (5-1) 
Hin 


It is called the wave impedance or the specific impedance and enters in the expression 
for the average power flow per unit area in the direction of the guide 


Ws = 3EvnHin = 3KnHinHen (5-2) 


The reciprocal of this impedance is the wave admittance 


M,, —_ ’ (5-3) 





the power flow is then 
We = MoE inEvn (5-4) 


In wave guides with perfectly conducting walls the various transmission modes 
carry power independently of each other. The field patterns are “orthogonal” to each 
other and may be “normalized”; that is, the transverse intensities for a typical mode 
may be expressed as follows 


E..n = V.F,(u, 2), Hin = I,F,(u, 2), I; = M.V., (5-5) 


f [F,.(u, v) |*dS 
ff Fa v)F,,(u, v)dS 


u and v are suitable coordinates in the transverse plane of the wave guide and the 
integration is extended over the entire cross-section. The coefficients V, and J, may 
be called respectively the normalized voltage and normalized magnetomotive force or 
normalized current associated with the mth mode. 

Calculating the total power carried in the mth mode, we obtain 


where 


i, 
(5-6) 


0, if mn, 


W = 1x «Ieee. (5-7) 


“ 


Thus, if we express our transmission formulas in terms of normalized voltages and 
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currents, the same impedance coefficient appears in the alternative expressions (3-1) 
for power and this impedance is also the ratio of the normalized voltage to the 
normalized current. 

Before going on let us see just what the above formulas mean in one or two special 
cases. Consider a wave guide consist- 
ing of two parallel metal strips of 
width a, separated by distance 38, 
hesetiaeaea Fig. 1. In the dominant mode the 
electric intensity is perpendicular to 
aia the metal plates and is distributed 

almost uniformly except near the 

edges and in the external region where 
the field is weak and little power is carried by the wave. Neglecting the edge effect, 
we shall assume that the electric intensity is constant 














Fic. 1. 


E, = Eo. (5-8) 


The normalized distribution pattern is given by 





1 
F(x, y) = seoeett (5-9) 
a/ab 
and, therefore, 
E, = VoF (x, ¥), Vo = Eqv/ab. (5-10) 


The wave impedance for transverse electromagnetic waves is Ko=~/u/eand therefore 


r 


ar 
Hy =IcF (x,y), Io = Hov/ab = ae (5-11) 


0 


In air Ko=approximately 377 ohms. The transverse voltage V between the plates 
and the longitudinal current J are 


V = Ey = Vov/b/a, I = ay = Iw/a/b; (5-12) 
consequently the characteristic impedance on the voltage-current basis is 
VebVvo 6b | 2 
Ky. = — = — = — Ko. (5-13) 
I alo a 


For the total power flow we have 


‘ +? we a? * — ‘ 
W = LK ol olo = 4M ob oVo = $Ky rl OE come (5-14) 
2Ky.1 
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so that in the present case = ~eaaiaelian a aaa 
b 
b 
Kw. = Kw = Ky = — Ko. (5-15) \ le A Ae ee 
a a t ~~ 
Consider now the 7E;,9-wave in a BY ------ 
rectangular wave guide, Fig. 2; for 
_ . ~ ne Fic. 2. 
this wave the field is given by 
~ a] . WX . Wx 
E, = FE, sin—» H, = H,sin—» 
a a 
(5-16) 
acs h2 \-1/2 
EF, = Kif,, Hy, = ME, Ky = /n/€ (1 a =) ° 
4a? 
The normalized field distribution function is 
2 wx 
F(x, y) = — sin —»> (5-17) 
ab a 
so that ) 
E, = ViF\(x, y), H, = IF \(x, y), 
cama niikienie (5-18) 
Vi = Ey/ab/2, I, = Hy/a/2b. 


In this case the maximum transverse voltage V across the guide and the total longi- 
tudinal current J are given by 


V = Eyb = Vy\/2b/a, 


, (5-19) 
I -{ Hdx = — Iy/ab/2. 
0 TT 
From these equations we have 
Vb 
Ky. =—=—K.. (S-20) 
I 2a 
The power transfer is 
* + * 
W = 3K; = 4MiViV;; (5-21) 
consequently , 
3b 2 
Kw, oe Ki, Kw,v ——— Ky. (5-22) 
8a a 


Now let us see what happens when we join two wave guides, each consisting of 
two parallel metal strips. Suppose that the frequency is so low that we do not have to 
worry about higher transmission modes. At the junction the transverse voltage and 
the longitudinal current must be continuous. This requirement is responsible for re- 
flection unless the characteristic impedances of the two guides are equal. The coefh- 
cients of reflection and transmission depend on the impedance ratio Ky,;/Ky,, of the 
two wave guides. As we shall find later the effect of the geometric discontinuity can 
be calculated equally well by concentrating attention on normalized voltages and 
currents. With respect to these variables the characteristic impedances of the above 
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wave guides are equal; but at the wave guide junction there will exist an effective 
ideal transformer with the impedance transformation ratio equal to Ky.;/Ky,;. In 
the case of ordinary low frequency transmission lines we prefer to think in terms of 
total voltages and currents; to think in terms of normalized voltages and currents 
would be to make simple matters complicated; but it will presently become evident 
that, in general, it is advantageous to introduce the normalized variables at least in 
certain stages of the analysis. 

Take an iris in a rectangular wave guide. We know that for frequencies between 
the lowest cut-off frequency and the next higher, the iris can be represented as a 
shunt susceptance. The value of this susceptance will depend on its definition; but 
the ratio to the corresponding characteristic admittance of the guide is an invariant. 
It is this ratio that appears in transmission formulas involving lumped elements in- 
serted in a uniform transmission line. If, however, the iris is between circular and 
rectangular wave guides, the ratio of the characteristic impedances of the two guides 
will also be involved and this ratio depends on whether both impedances are defined 
on the power-voltage basis or the power-current basis. It is evident, therefore, that 
in this case the iris cannot behave as a simple shunt susceptance. The theory which 
we are now evolving permits us to prove that in the more general case the equiva- 
lent transducer for the iris consists of two shunt susceptances, corresponding to the 
two faces of the iris, and an ideal transformer between them. The transformer ratio 
depends on the particular voltage-current set we happen to choose for our analytical 
work but our final transmission formulas will be independent of this choice. The 
degree of arbitrariness involved in the choice of “V” and “J” is of the same kind as 
that involved in the choice of coordinate systems or of units. In elementary analysis, 
a particular choice was so natural that a mistaken notion spread abroad that this 
choice was a necessary one. 

6. Aniris between two wave guides. Let us now obtain an exact equivalent circuit 
tor an infinitely thin perfectly conducting iris between two wave guides of arbitrary 

















“(1) —_ i 
(2) 
| —_ — > 
- 77 | ; 
Fic. 3. 


cross-section (Fig. 3). The constants of this circuit depend on the particular trans- 
mission mode under consideration; that is, there is one equivalent circuit for transi- 
tion from each transmission mode in wave guide 1 to each mode in wave guide 2. The 
most important case is that of transition from the dominant mode in one wave guide 
to the dominant mode in the other, and in the following analysis we shall keep this 
case specifically in mind; but the analysis applies to any other case. We shall use 
Cartesian coordinates in our equations; but this does not mean that our analysisis 
restricted to rectangular guides. 

Suppose that the transverse field of the incident wave at the surface of the iris is 
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Ei(x, y) = ViF1(x, y), 
Hi(x, y) = MiViF (x, y), 


(6-1) 


where V} is the normalized incident voltage. In response to this impressed field, we 
shall have some field over the aperture of the iris. Let f(x, y) be the tangential electric 
intensity over the aperture; then in wave guide 2 the “transmitted” tangential elec- 
tric intensity is 
Ei(x, y) = f(x, y) over the aperture, 
(6-2) 
= (0) over the screen. 
In wave guide 1 the total tangential electric intensity, that is, the sum of the incident 
and the reflected intensity, must be 
ViF (x, vy) + Et(x, vy) = f(x, y) over the aperture, (6-3) 
. 6-3 
= 0 over the screen. 
The function defined by (6-2) may be expanded into a series of normalized or- 
thogonal functions appropriate to wave guide 2; thus 


Ei(x, y) = >> VF, (x, y). (6-4) 


n=1 


The tangential magnetic intensity is then 
Hii(x, y) = > M.V,F,(x, 9). (6-5) 
n=1 


The function defined by (6-3) can be expanded into a series of normalized orthogo- 
nal functions appropriate to the wave guide 1; thus 


ViF (x, y) + E(x, y) = > Vil n(x, y). (6-6) 


n=1 


The reflected tangential intensity is therefore 


Ei(x, y) = (Vi — VDF A(x, y) + DS VF, (2, 9). (6-7) 


n=2 


The corresponding tangential magnetic intensity is then 


Hix, y) = — MAVi — VDFA«, ») — SS MeV nF (2, 9). (6-8) 


n=2 
The transfer of complex power through the aperture must be continuous; therefore 


rity . r tr* — oe f 
2M,ViV: — > MV .Vs = >, UV. (6-9) 
n=1 n=1 
The voltage reflection coefficient gy is defined as the ratio of the reflected voltage 
V,— V; to the incident voltage V{; it may be obtained from (6-9) if we divide the equa- 


tion by 2MiViV¥#; thus 
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° * ~~ a =a* ~ 
1 vi “ ViVn M,\ViV o M.V,V, 
- as [14 2 = |+ Vis 1 re ar: | (6-10) 
2 1 








1 + qv V; n=2 M;, VW 2M,V iV} n=2 MV, 
Consider now the complex power flow into the second wave guide 
*e logs lagcso 
WwW = 4 MV, + — >> M.V.V, - (6-11) 


The form of this expression is such that from the input end various transmission 
modes appear to be in parallel. It is not exactly that the characteristic admittances 
Mi, Mz, Ms, --- are directly in parallel ; if we select the first admittance for reference, 
the others are transformed in the ratio V,V,*/ Vi V# before being connected in parallel. 
In any case the net effect on the input admittance is the same as would be obtained 
if we had an admittance Y in shunt with a transmission line maintaining only the 
dominant mode. Thus we can write 





W = 430,,V,V5+ 490,77, (6-12) 
where 
ns t> T7* 
Y= > uv,——- (6-13) 
p> ViVi 
The ratio of the shunt admittance to the characteristic admittance 
J a o M00. 
ei ee (6-14) 
MM, << M.V,V 


is an invariant. It has the same value regardless of a particular basis for definition 
of admittances and it depends only on the form of distribution of the tangential 


electric intensity over the aperture. 
Similarly for the admittance ratio looking from the iris into wave guide 1 we have 





MW Ve 
ag MiWiV* 


(6-15) 


iM 


If the frequency is in the interval between the lowest cutoff frequency and the 
next higher, then M2, M3, - ++ arereactive and the shunt admittances are pure sus- 


ceptances 
Y=iB, Y= iB. (6-16) 


For frequencies higher than the second cutoff frequency an iris entails some power 
loss to the dominant wave. The lost power is contributed to one or more higher 
transmission modes. This is analogous to what happens when a doublet antenna is 
inserted in shunt with a parallel pair or at the end of it. The plane wave guided by 
the parallel pair loses power; this power is then carried away by a spherical wave 
which originates at the junction. One mode of energy transmission is partly trans- 
formed into another. Usually there is also an energy exchange between a local field 
and the plane wave; this results in reactance. At lower frequencies an ordinary coil 
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(or a capacitor) inserted in a transmission line acts just like an iris; electrically it does 
not matter just what physical means we happen to provide for a local storage of 





energy. 
We now can rewrite (6-10) as follows 
pits ‘ia (1 we =) t i nn* (1 4- =) (6-17) 
1 + qv 2 M, 2 M, M, 
where 
nn* = Lie (6-18) 
ViVt 


The reciprocal of the voltage transmission coefficient is 
1 Vi 1 vs; 171 Yy 1 M, Y 
— = SS = —] — [14+ —) + — ert* — (14+) ]. (6-19) 
pv V; 1+ qv Vi n 2 M, 2 M, M, 
It is a simple matter to prove? that for transmission lines coupled as indicated in 
Fig. 4, py and gy are given precisely by equations (6-17) and (6-19). In Fig. 4 the 


i:n* 





“~ 


M 





Fic. 4. 


transformer ratio 1: n* is indicated for the impedances rather than for the admittances 
in order to conform to the established practice. If m=1, which is always the case when 
the wave guides on both sides of the iris are the same, the admittances Y and Y of the 
two faces of the iris are just in parallel, and the transformer can be omitted. 

The exact numerical values of n, ¥/M:, Y/M: are found by solving the appropri- 
ate boundary value problems.* The approximate values can be obtained quite easily 
if we assume a reasonable form of distribution of the tangential electric intensity 
over the aperture,‘ and of course, we can always calculate these quantities from 
measurements of the transmission and reflection coefficients for waves moving from 
one wave guide into the other. Thus, 


os 


l+o7 M1 1+ q7 Ki 
ial ees tae ie cone ey (6-20) 
i + gy M, i+ Ty Ky 
where gy is the voltage reflection coefficient for a wave moving from left to right 
and gj is that for a wave moving in the opposite direction. 


2 See for instance S. A. Schelkunoff, Electromagnetic waves, D. Van Nostrand Company, Inc., New 
York, 1943, p. 212. 

* For example, see S. A. Schelkunoff, The impedance of a transverse wire in a rectangular wave guide, 
Quarterly of Applied Mathematics, 1, 78-85 (1943). 

4S. A. Schelkunoff, Electromagnetic waves, p. 491. 
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If the iris is not indefinitely thin, there is a section of a wave guide between the 
two faces of the iris. 

While the iris acts effectively as a lumped impedance, the field associated with 
it is actually distributed. Even if the frequency is such that the iris is reactive, the 
field extends to some distance on either side of it. Near the cutoff for the second trans- 
mission mode this distance may be quite large; but ordinarily the field extends roughly 
to a distance comparable to the transverse dimensions of the guide. There will exist, 
therefore, a mutual impedance between those faces of two nearby irises which face 
each other. For frequencies above the second cutoff, the mutual impedance may, and 
usually will, exist even between two distant irises. All these considerations do not 
affect our essential picture of electrical properties of wave guide discontinuities; they 
affect merely the numerical values of various impedance and admittance functions. 

In the above equations we have treated E, and H, as if they were scalars; in gen- 
eral, they are vectors. However, the analysis is similar to the above and the final 
formulae are the same. 

In the case of coaxial pairs or wave guides formed by parallel metal strips the 
dominant wave is transverse electromagnetic. If the edges of the iris are normal to 
the lines of force for the dominant wave, the voltage between the edges is equal to 
the transverse voltage across either guide; the total voltages associated with higher 
transmission modes are equal to zero; and the transformer ratio is unity provided we 
base our transmission diagram not on the normalized characteristic impedance but on 
the conventional impedance K which in this case equals Ky,1, Kw,1 and Kyw,y. 

7. Reactances in series with wave guides. An example of a reactance effectively 
in series with the wave guide is shown in Fig. 5 which represents a circular wave guide 








Sat t 
Stat — 2a 
Pe { 
It 
Fic. 5. 


and a narrow radial transmission line.’ Let us suppose that we are concerned with 
transmission of a 7 Mo,:-wave. For this wave the field is circularly symmetric. Mag- 
netic lines are circles coaxial with the tube, and electric lines are in radial planes. 
It is practically self-evident that the radial line is in series with the guide, and that 
in parallel with the radial line there is an impedance associated with the gap. If the 
frequency is between the lowest cutoff frequency and the next higher, this “gap im- 
pedance” or fringing impedance is capacitive and is 
of little importance except when the impedance of 
the radial transmission line is high. For frequencies 





Ke above the second cutoff, the gap impedance is in 

z part resistive on account of power transfer from the 
Ks dominant wave to the higher order waves. As seen 

from the gap, the impedances of various waves in 





the guide and the impedance of the radial wave are 
in parallel; the two halves of the guide are in series; 
Fic. 6. and the impedance diagram looks like that shown 





5S. A. Schelkunoff, U. S. Patent 2,155,508, April 25, 1939. 
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in Fig. 6. The same diagram is shown in Fig. 7 where the characteristic impedances 
K, and K, have been “expanded” into semi-infinite transmission lines; the impedance 











consisting of 2Ke, 2Ks3, - - - in parallel is represented simply as the gap impedance Z,. 
2 
Fic, 7. 


Starting with equation (10.17-1) of “Electromagnetic Waves,” we can obtain the 


approximate gap impedance by the method explained there. In this case, however, 
the following elementary derivation is preferable. To begin with, let us remove the 
radial line and assume that the electric charge is being transferred across the gap by 
an impressed voltage V‘. The total conduction current J in the tube is the sum of 
currents associated with the various transmission modes. Thus for the input current 
we have 


T=ht+I+Is+---. (7-1) 


The input power is then 


W = 3Vi" = BVT + AVID+--- 
7-2 
= $i + iehele t ++: , — 
where Z;, Z2, +++ are the input impedances of individual waves; that is, 
Limes yt (7-3) 
Fe tae nen Fa may coe. Das 
_-_ 


In the above equations we have tacitly assumed that the gap is very small and 
the current associated with each mode does not vary in the gap. This restriction will 
presently be removed. The total power contributed to the wave is divided between 
different modes; one-half of it is carried to the left and the other half to the right. 
The power carried in one direction in the mth mode is }K{r/,J,*; thus we have 


Z, = 2p. (7-4) 


Actually the applied voltage is distributed in the interval (—s/2, s/2) around the 
midpoint z=0. Assuming that the distribution is uniform, we may write the con- 
tribution to the total current associated with the mth wave at point z due to an 
elementary voltage at point 2 as follows 


A 


,, a2 
T,eTale-#l —, (7-5) 
: 
where J, is the amplitude at the source. The total current at point z is then 


A 


F. s/2 ; dé 
ae | ee (7-6) 


5 3/2 s 
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The power contributed to the mth wave is then 


a/2 dz 
Wa vi f I (s —= (7-7) 
—s/2 S 
Thus we shall have 
Wn = 4xnViTny (7-8) 


where (assuming that I, is real) 


1 pet pen 2p 2 Mtr 
mS dz eln!? 7ldz = = = eS 7-9 
X 5 ’ 19.9 

os —s/2 —s/2 I nS Iss? 


On the other hand 


. ,(n) t 
VW » = Ks tT: (7-10) 
consequently 
yi 2Kwr ‘ 
——— : (7-11) 
Pn Xn 


Since x, decreases with increasing n, the successive components of the gap admittance 


; a ae ‘ X2 X3 X4 a 
re eo Fee Fee 28 Se KB ARC (7-12) 
2) 9 k (3) 9 kA) 
2Ky, 2Ky, 2Ky, 
decrease. 
The typical K, is given in problem 8.10 on page 509 of “Electromagnetic Waves” 
ia P ; k? Ag? "1/2 - 
aw 27 = . (7-13) 
4riwe a? ? 
where &, is the mth zero of Jo(x). For sufficiently large n, therefore, we have 
2 _(n) r,s k,,s 
Z, = ',Kwus = ie + (7-14) 
4triwe 4ariwea’* 
The impedance of the radial line is approximately 
F 2rl 
Z = 60: tan —- (7-15) 
a r 


A more accurate expression in terms of Bessel functions may be found on page 269 
of “Electromagnetic Waves.” 

8. Conclusion. The ideas developed in this paper are adequate for expressing 
transmission properties of wave guides with discontinuities in terms of impedances 
and admittances associated with these discontinuities. These impedances are reactive 
if the frequency is such that the energy in either guide can be transmitted to any 
distance in only one mode; otherwise, the discontinuities present some resistance for 
the mode under consideration and a negative resistance to those other modes which 
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participate in transmission of energy. The finding of exact values of impedances re- 
quires solution of corresponding boundary value problems; but frequently good ap- 
proximations can be found by making reasonable a priori assumptions on physical 
grounds. In fact, the point of view outlined in this paper makes it easy to make such 
assumptions. 

More complex discontinuities can be analyzed into simpler discontinuities. The 
discontinuity shown in Fig. 8 is equivalent to an ideal transformer between two wave 
guides; across the left “winding” of which there is a small shunt capacitance® and 
across the right winding there is the capacitance’ associated with the annular disc 


eae 











Fic. 8. 








Fic. 9. 


looking into the second guide. In parallel with the latter capacitance there is the series 
combination of the impedance of the radial line and the second guide itself. We may 
express these ideas by the diagram shown in Fig. 9, where the inductance is used to 
designate the radial transmission line only because this line, when it is short, is ap- 
proximately an inductance. 

More generally, the discontinuities should be represented by impedances distrib- 
uted along the guide, as in fact they are. Finally, the section of the guide with the 
discontinuities may be replaced by an appropriate 7-network. 

Recently J. R. Whinnery and H. W. Jamieson® have obtained explicit expressions 
for the capacitances of numerous types of “step discontinuities” in transmission lines 
formed by parallel conducting planes. They show how to apply these results to coaxial 
conductors. They find theoretical predictions in good agreement with measured 
values. The equivalent circuits given by Whinnery and Jamieson do not contain ideal 
transformers; this is because for transmission lines comprised of two conductors, the 
transformer ratios at discontinuities are equal to unity and the transformers may be 


omitted. 


6 In the first approximation this capacitance may be neglected. 
7 We assume that we are operating below the second frequency cutoff; otherwise there will also be a 


conductance. 
8 J. R. Whinnery and H. W. Jamieson, Equivalent circuits for discontinuities in transmission lines, 


I.R.E. Proc., February 1944, pp. 98-114. 








THE DISTORTION OF THE BOUSSINESQ FIELD DUE 
TO A CIRCULAR HOLE* 


BY 
A. BARJANSKY! 
The Brush Development Company 


1. Introduction. One of the most important problems in the theory of elasticity 
is the solution of the biharmonic equation Y°¢ =0, where ¢ is Airy’s stress function, 
for a given group of boundary conditions. As is well known, the most common ap- 
proach to the solution of this problem consists in selecting a system of coordinates 
particularly suited to the region studied. 

Thus, using bipolar coordinates, G. B. Jeffery has given the general solution of 
the plane problem, that is, of the biharmonic equation in two dimensions, for regions 
bounded by non-concentric circles (Ref. 1). A clear, but not quite complete, treatment 
of Jeffery’s method can be found in Coker and Filon (Ref. 2). This method has re- 
cently been used by R. D. Mindlin for the determination of dead loads on tunnels 
(Ref. 3). 

The present paper is an attempt to apply Jeffery’s approach to the problem 
of the distortion introduced in the so-called plane Boussinesq field by the presence 
of a circular hole. Starting with the stress function ¢ of the undistorted Boussinesq 
field, an auxiliary stress function x will be found such that =@+y satisfies the differ- 
ential equation and all the boundary conditions. The stresses and strains in the dis- 
continuous field can then be directly determined from the derivatives of ®. 

2. The Boussinesgq field. Boussinesq and Flamant have given the solution of the 
biharmonic equation for the case of an isolated force P acting at a point on the bound- 
ary of a semi-infinite plane. Their solution, which can be found in all standard texts 


(see, for instance, Ref. 4, p. 82) is: 


P 
o; = — —76 sin 0 (1a) 
T 
for the case of a normal force, and 
| 
od: = — —récosd (1b) 
T 


for a force parallel to the boundary. The significance of the symbols is shown in 
Figs. 1a and 1b. 

In the simple Boussinesq problem, the only boundary conditions are that the 
stresses, both normal and shearing, must vanish along the straight boundary (except, 
of course, at the point of application of the force) and also must tend to zero as one 


* Received Sept. 23, 1943. 

1 The writer wishes to express his thanks to the following persons who have assisted him in the 
preparation of this paper: Professor M. S. Ketchum of Case School of Applied Science, who suggested the 
problem; Dr. W. N. Dudley, also of Case, for several helpful hints; and especially Dr. H. G. Baerwald 
whose assistance in many mathematical details has been invaluable. 
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moves away from the point of application within the half-plane. With the appearance 
of the circular discontinuity (Fig. 2), the above conditions remain, and a new one is 


Pe p 


bh t . 








Fic. la. Fic. 1b. 


added, determined by the nature of the discontinuity. Thus, if it is a hole, both normal 
and shearing stresses must vanish along its periphery. 








The region is thus bounded by a 
circle and a straight line; the latter can * 1p 
be considered as a circle of infinite y 


radius, so that here is a case of a region 
bounded by two non-concentric circles, ag 
to which Jeffery’s method is applicable. ¢ 
3. Bipolar coordinates. Jeffery’s : r 
method consists essentially in intro- 
ducing a system of curvilinear coordi- 
nates, called bipolar coordinates in 
works on elasticity. Two poles, A and 
B (Fig. 3) are taken at abscissas +a Fic. 2. 
along the X-axis, and the location of 
any point is determined with respect to these poles by the quantities 








rT) 
& = log — n = 0, — 62. 
2 
The lines »=constant are circles passing through A and B, while £=constant are 
a system of circles with centers on the X-axis. Some of these lines are drawn on Fig. 3. 
If a circle of diameter d has its center h units from the horizontal axis (Fig. 2), 
it is easy to show (see Ref. 1) that the proper polar distance a is determined from 
a? =h?—d?/4 and that the value £ of £ corresponding to the circle is §)=cosh— 2h/d. 
The cartesian coordinates can be expressed as follows in terms of the bipolar: 
a sin 
a sinh & sin (2) 


’ y . 
cosh § — cos 7 cosh § — cos 7 





When the biharmonic equation is expressed in bipolar coordinates, it is found 
convenient to write it, not in terms of the usual stress function x, but in terms of 


x/J, where J has the value 
a 





cosh § — cos 9 


the stresses are also expressed as derivatives of x/J. 
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The bipolar solution for x/J used by Jeffery has the general form: 
x/J = Bt cosh § + (— BE +G cosh 2¢ + H sinh 2 + F) cos n 
+ (G’ cosh 2§ + H’ sinh 2§ + F’) sin n 


+> { [E, cosh (& + 1)€ + Fy sinh (& + 1)E + Gy cosh (k — 1)E 
k=2 


+ H;, sinh (k — 1)&] cos kn + [Ed cosh (Rk + 1)E + Fi sinh (& + 1)E 
+ Gi cosh (k — 1) + AZ sinh (k — 1)é] sin kn}, (3) 


where all the B’s, E’s, F’s, G’s and H's are constants. This series will be assumed 
convergent and differentiable for the time being. 

Here the terms independent of 7 and those containing cos 7 or sin 7 are used ex- 
actly as they appear in Ref. 2 (Eq. 4.066 and paragraph 4.07), but those containing 
functions of multiples of 7 come directly from Ref. 1 (Eq. 21), with some slight 


changes in nomenclature. 

















” 


Y 





xv 
Fic. 3. 


4. General procedure. As was said before, the presence of the circular discon- 
tinuity causes a modification of the Boussinesq functions ¢: and ¢2 into ®,, and 
the latter having to satisfy the biharmonic equation and all boundary conditions 
Also ®;=¢1+ x1, P2=¢2+x2 where x1 and x2 are auxiliary stress functions of the gen- 
eral form (3). Now since both @’s and both y’s satisfy the biharmonic equation, which 
is linear, so do ®; and ®,. As to the boundary conditions, ¢: and @¢2 satisfy them along 
the straight boundary and for remote points. Therefore, x: and x2 must be so selected 
that: 
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(1) they give vanishing stresses for remote points (0, »-0); 

(2) they give zero normal and shearing stresses along the straight boundary; 

(3) in combination with the known functions ¢; and ¢: they satisfy the boundary 

conditions at the circular discontinuity. 

In the next paragraphs, conditions (1) and (2) will be considered first and their 
application will determine some of the hitherto arbitrary constants of Eq. (3). 
Then the function x satisfying conditions (1) and (2) will be added to ¢: or ¢2 (ac- 
cording to whether a normal or a tangential load is studied), yielding 


=¢+x, P. = oo + x. 


Finally the remaining constants of x will be determined in each case by the conditions 


at the inner boundary. 
5. First and second boundary conditions. The stresses are expressed as follows in 


terms of bipolar coordinates: 


2 


te] te) 0 x 
ao; = | (cosh £ — cos n) —— — sinh  — — sin » — + ccsh :|(*). 
On? 0é On J 


ae  @ . 4 x 
ac, = | (cosh £ — cos n) —— — sinh § — — sinn — + cosy |{—}, (4) 
ag? at an J 
0? x 
art, = — (cosh  — cos 7) ——{— }. 
0&0n \ J 
The first condition necessitates 
o: = 0, o, = 0, re, = O for (€, n) > 0; 


and the second 
o: = 0, te, = O for & = 0. 


The first condition is seen from Eqs. (4) to be equivalent to 


x 
— = 0 for (t,) 30 
P 


from which, immediately 
G + PF = 0, G = —F and E; + G; => 0, Gi == Ex. 
For the second condition, 7;,=0 for £=0, which yields 


B k+1 k+1 
H=—, AW=0, Mm=-- F,, Hi =- Fy; 
y) k—1 k—1 








and from o;=0 for §=0, 
Gi = = RY . 


Thus the stress function satisfying boundary conditions (1) and (2) assumes the form: 


= Bé cosh § — [B(t — sinh £ cosh £) + 2F sinh? £] cos 7 + G’ cosh 2¢ + F’) sin 9 


Ss 
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+ 2.5 7 { [E,(& — 1) sinh & sinh hE 


+ F,(k sinh — cosh kE — cosh & sinh E)] cos ky 
+ [E/(k — 1) sinh é sinh kt + Ff (& sinh £ cosh RE — cosh £ sinh é) | sin kn}. (5) 


6. Third boundary condition. The value of x//J from (5) is now added to ¢; (for 
normal load) or ¢2 (tangential load), and the remaining arbitrary constants of (5) 
determined by the conditions at the boundary of the hole, which are that both normal 
and shear stresses vanish on the periphery, i.e. (o¢);,=0, and (7¢_):,=0. 

In order to coordinate the functions ¢1, ¢2, on one hand, and x1, x2, on the other, 
the system of axes shall be selected so that the y=0 (or 7=0) axis passes through the 
center of the hole, as shown in Figs. 2 and 3. Then the concentrated force, whether 
normal or tangential, will act at a point y=vyo, and the stress functions ¢: and ¢2 
become 

F y— Yo P y— Yo 


od: = — — (y — yo) tan—' ———» go: = — — tan“!- 
T x 7 x 





Transforming this into bipolar coordinates (Eq. 2) one has 





F yo(cosh £ — cos n) — asin 17 yo(cosh £ — cos 7) — asin 1” 
¢, = — — ————  tan-? s——_____________ , 
T cosh F — cos 7 a sinh & 
P a sinh é yo(cosh ——cosn) — asinyg 
¢:=— — — tan-! ———_-____________ 
T cosh ¢ — cos 7 a sinh 


But, as was said before, in treating problems involving bipolar coordinates, it is easier 
to express stresses not in terms of the stress function itself, but in terms of the stress 
function divided by the quantity J, so that: 


yo(cosh £—cosn) —asinyn 





¢:/J = — [yo(cosh — — cos n) — a sin 9] tan~! ———————_ ; 
Ta a sinh 
r . yo(cosh  — cos n) — asin us 
¢:/J = — sinh & tan~) ———_—_—_______—_ 
T a sinh & 


These two expressions must now be written in Fourier series in 7 to be comparable 
with the auxiliary functions x: and x2 of (5). The coefficients of these series are found 
by means of the usual integrations which are presented in detail in the Appendix. 
The results are as follows: 


¢:/J = + 2 T; cos kn + U; sin kn), 


with 


te 
= ll B cosh e| tan (tan B coth &) + (= _ 8) sinh | — ef ——__- 
T 
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P 
T,; = — tan al « sin 6 cosh — — tan~! (tan B coth &) — (= _ é) sinh |, 
Tv 
r j tan B T ‘ 
1= ) E + (1 + e*) cos e| — tan~' (tan B coth &) — (= — 2) sinh : ; 
WT \ 2 


and for k=2 


T,.=—- win (—1)* wee dei avr. sinh £—cosh *f1+(— cos =|} 











r | k?—1 cos B 
Pett tan B cos k k sinh £—cosh tan B sin kB 
jam PO RO SOR a, 
r | d k?—1 k cos B 
Ro k=o : 
$:/J = — + > (Ri cos kn + S; sin kn), 
é k=] 
where 


Ro FPsinhé T F 
-_ = | tan- (tan 6 coth &) + (= _ 8) sinh ‘|, 
2 7 2 
Psinhé e* ’ 
R, = —- = (— 1)* sin &£, 


Tv 





P sinhé e* 
S, = ———_ — [1 — (— 1)* cos A]. 
vg k 
In these expressions B=tan-! yo/a (Fig. 2). The above formulas refer to the case 
80 (see Appendix). For the important case B=0, i.e., yo=0 


To P 
oc". ? oe ef, T; aa 0, Ui; ag 0, 
2 T 
for k=2 
2P e**(k sinh — — cosh &) 
T,=- —_—_——— —— for k even, T. = Ofor k odd, U; = 0, 
T k? —1 
and 
Ro 2P e*é ’ 
a 0, R, = 0, Ss, =— = sinh £ for k odd, S; = 0 for k even. 
Z TT 


The following special cases will be considered in the next section on applica- 
tions: (A) half-space containing a hole bounded by £=£ and subjected to normal 
load; (B) same region subjected to tangential load. 


APPLICATIONS 


A. Hole subject to normal load. Here 


®,/J = o:)/J + x:/J, 


or 
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1+cos =| 
2 cos B 


D 


®,/J = Bé cosh £+— jtan B cosh e| tan (tan B coth£)+ (=-s)sinn | — e&— 
T v4 





~ 4 5( (¢—sinh ~ cosh £)+2F sinh? é 
tc T ; \ 
——tanB E sin B cosh €— tan! (tan 6 coth £)— ——- sinh £ f cos B 


T 
P [tan B 
+ \e cosh ete +—| =F { 1+(1+ 7) cos B} 
T 
— tan! (tan B coth »-(<-8) sinh | sin B 


“ 


ak 


1 
nw 
a 


nr (zal —1] sinh é sinh kE+F;,,[k sinh cosh k§E—cosh € sinh hE | 


k 











k tan B sin k k inh —cosh cos kB 
ans )* ee ae E = “= [14+(-1) |b) cos kn 
2 R-1 cos B f 
- 2 
+ i fe yaa ( 2 [fe | sinh & sinh AE+F{ [& sinh — cosh kE—cosh & sinh £| 
k=2 








k?—1 k cos B f 


The condition (7;,):,= —(cosh —cos 7)d?/d£0n(@/J)=0 amounts to equating 
to zero at £=£, the derivative with respect to — of each term except the one inde- 
pendent of 7. As to (o;);,=0, this can be shown to require that for k=2 each term 
be zero at £=£. Thus, for each term, two equations are available; and this is sufh- 
cient to find all of the remaining constants, with the exception of F’ in the term in 
sin 8. The constant F’ remains indeterminate, and can therefore be taken as equal 
to zero. By solving the two equations for each term, the following values are found 
for the constants: 


P 2 cosh? & 2 sin? B(cosh? 95+ 
B= : {ts tan 8 sin p(t) aK +e 
2m sinh? & sinh £5 cosh? £)—cos? B 
1 + cos B\ 


—(mx—28) tan B cosh (cosh? §:—}) —coth & —————;7,, 
cos B f 


k—1 é k k sinh &—cosh £ [tz in kB}) 
-—— Pett} (- 1)* oe meosh gan | (ay =) sin kn. 


2a 











je tan B sin B cosh £)—tan B sinh £9 cosh £9 


sin 8 cos B wT 1+ cos B 
cosh? £,—cos? B 2 2 cos B 


—P ; sin 8 cos B T 
G’ =—————| er sin B+ —_________— (=-s) cosh £ |, 
2m sinh 2£ cosh? §)—cos? B 2 


F’=0, 


~ 26 sinh? Eo 





P 
Ey = — = {(—1)* tan B sin kB(k sinh £) cosh £+sinh? &)—e*te sinh kéo| 
Tv 
+ [1 — (— 1)* cos kB/cos B]k sinh? £0} [sinh? kEo— &? sinh? >|, 
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P 
FF, = — I (—1)*(k—1)k" tan B sin kB[k sinh & cosh f)—e**s sinh kéo| 
T 


+ (k+ 1) iff +( — 1)* cos kB/cos B|[R? sinh? £) — & sinh £ cosh &) — e~**» sinh kto|} 
-[sinh? kEo — &? sinh? |", 


? 


Ei =— (—1)*tan B cos kel k sinh k§) cosh £)+sinh? §)—e*fe sinh hEo| 


2r 
—[k- tan B + (—1)* sin kB/cos B]k sinh? £o]} [sinh® kEo — &? sinh? fo], 
P 
Fj = : -}(—1)*(k—1) Rk“ tan Bcos kB [k sinh £9 cosh £9+e** sinh kEo| 
PAs 


~(k+1)—[k- tan B+(—1)* sin &8/cos B|[k? sinh? {)— & sinh & cosh £o—e* sinh kEo} } 
-|sinh? k&o — # sinh? &)|-. 


To test the suitability of this expansion, it is sufficient to examine the terms of 
the auxiliary functions (Eq. 5) for the values of the constants given above. The co- 
efficient of cos kn in the general term of the latter equation is seen to consist of two 
parts, one multiplied by (—1)* and the other not. The first part forms an alternating 
series the general term of which tends to zero, so that the alternating series is con- 
vergent by a well-known theorem. The second part is found to converge outside the 
circle £=£ by the ratio test. The same is true for the coefficient of sin kn. Thus 
the above expression for ®,/J is a uniformly convergent series in 7 in the region con- 
sidered. 

Because of the great complexity of the expression involved, only the case B=0 
will be considered in more detail. For that case 

P cosh & P 
Bo = — — ———- Fo -— —_ —e 


' ‘ ~ a, G’ = KF’ = 0, 
x sinh® & 2m sinh? & 








. P k sinh? £5 ™ 
Ey’ = -— — ————"—_— for b even, E;” = 0 for k odd, 
mw sinh? kip — k? sinh? £o 
P k*sinh? £ — k&sinh £9 cosh &) — e**o sinh k 
FO = — Paacsnedl i > ec ten, Ot « Gi bell, 





. (k + 1)(sinh? ky — &? sinh? £) 
EX = Fi = 0. 


The stress function becomes 





“ P cosh & P 
@f /J=- - — cosh ¢ — -—e§ 
7 sinh? & T 
P T cosh & sinh? & 
+. ——— (¢ — sinh £ cosh £) + — | COS 9 
T sinh® £ sinh? £o 





2P k(k — 1) sinh? 
| \ = sinh é sinh ké 


pag ™(Rk — 1) Lsinh? kt — k? sinh? & 
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k? sinh? £) — &sinh £) cosh &) — e*€e sinh k&o ; . 
- —______—____4!— (ksinhé cosh kt — cosh £ sinh ké) 


(k + 1)(sinh? kkp — &? sinh? &o) 
ek 
+ —— (ksinh ¢ — cosh £) | cos kn 
k+1 
with the summation extending over even values of k only. 
The most significant stress is the hoop stress ¢, at the periphery of the hole, = £o. 
Substituting the above value of #;/J into the second of Eqs. (4), the following series 


is obtained: 













P P 2P cosh & ' ae ; 
ac, = — — (1 + coth? &) — —- ———— cosy + — (5 coth? && — 1) cos 2y 
7 mw sinh? & T 
te 2 sinh £ sinh 2£ 4 sinh & sinh 4£ 
oo” __—_—_— + —___— ———— | cos 3n 
x L sinh? 2&) — 4 sinh? & sinh? 4&5 — 16 sinh? & 
4P 4 sinh £ sinh 4£ 
ee Oe 9 
x sinh? 4& — 16 sinh? £o 
ZF 4 sinh £ sinh 4&5 6 sinh & sinh 6£€ 4 
mabe ee meme ee gee ** 
x Lsinh? 4&) — 16sinh?& ~— sinh? 6) — 36 sinh? & 
40, 
| - za 
ry cosh ¢ = &h 3 
d 
5 | 
R 
+") 
€ 2 
Rn | LIMIT CURVE 
a] F-+ 0 
v } 
S J 
i 
80 100 
-2 ANGLE & IN DEGREES 
4| 
-6 





Fic. 4. 


Fig. 4 is a graphical illustration of the above formula. In that figure, the “stress 
factor” is plotted for different values of 2h/d=cosh &. By “stress factor” is meant 
the ratio of the stress a, to the stress which would have existed under the same load- 
ing at a point corresponding to the center of the hole, if the latter had not been 
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drilled. If there had been no hole, the point corresponding to its center would have 
been under a stress —2P/mh (compression) so that the stress factor is the ratio 
—go,/(2P/mh). Therefore, a positive value of the stress factor represents compression, 
a negative value, tension. 

It is seen from the figure that for each curve there exists a tension directly under 
the load (a=0), which becomes a compression as a is increased, reaches a maximum, 
then decreases, and becomes tension again when a approaches 180°. For low values 
of cosh £, i.e., of the depth-to-diameter ratio, there exists a secondary maximum 
of tension in the neighborhood of a= 20°. 

As cosh £ increases (as the hole gets deeper and deeper), the stress factor 
curves tend towards the “limit curve,” which is simply the graph of 1—2 cos 2a. 
The latter expression (Ref. 4, p. 77, second of Eqs. (58)) is obtained by assuming the 
hole to be in a field of uniform compression, equal to the compression —2P/zh at the 


center of the hole. 
B. Hole subjected to tangential load. Now the total stress function has the form 


2/J = g2/J + x2/J, 


where ®, is the total stress function, ¢2 is given by (1b) and xz is of the general 
form (5). The heretofore arbitrary constants are determined by the conditions at the 
inner boundary, which are the same as in the preceding case. The remaining constants 


are found to be: 

















P sin B sin 2B coth £9 cosh? & 
B= |= —— + ——___—__ — (sz — 28) ——— |, 
2m L sinh? £5 cosh? £5 — cos? B sinh & 
P Te®esin sin B cos 
as fc JA css (= m 2) —_ to 
2rL sinh & cosh? £5 — cos? B 2 
G’ P e7ko (1 4 8) 
,= = — — cos B), 
2m sinh 2£ 
F’ = 0, 
R P (— 1) sin k8 k? sinh? &) + & sinh £ cosh &) — e*€e sinh kéo 
sia 2x k sinh? kt) — k? sinh? £o 
P : sinh? & 
F, = — —(— 1)*(k — 1) sin kp : : ’ 
2r sinh? kf) — k? sinh? £5 
or P [1 — (— 1)* cos kB| sinh? &) + &sinh £ cosh &) — e*e sinh kéo 
De k sinh? k&y — k? sinh? & ; 
sinh? £ 





P 

7!) = fe im ows k ° 
we 2r Cf at nee #8] sinh? kip — k? sinh? £o 

The resulting Fourier series can be shown to converge as in the previous case. 

C. Conclusion. In the above paragraphs, a method was presented for computing 
the distortion of the original Boussinesq field when a hole is introduced. Other inter- 
esting results can be derived by simple means; thus, by superposing on the above 
stress functions ®; or ®; one of the solutions presented in Refs. 1 and 2, it is possible 
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to obtain the stress system for a Boussinesq field containing a hole, the periphery of 
which is subjected to a uniform pressure. Another extension of the above method, 
on which the writer is working at present, can be used to solve the case of a Boussinesq 


field containing a rigid disc. 


APPENDIX 
The decomposition of the Boussinesq stress function into a Fourier series in 7. 


We shall begin by decomposing the shear stress function 


2 P- yo(cosh — — cos n) — a sing 
— = — sinh ¢ tan~! —————__ —_———- 
T a sinh 


The different Fourier coefficients are given by 


} 1 2r do 
Ry + 1S; = f — e**dn. 
ae 


T 


This can be simplified by introducing the angle 8 = tan! yo/a. 


2 P yo(cosh £ — cos) — asiny 
— = — sinh ¢ tan-) ——_—__—_______—__—_- 
e T a sinh & 
a sin 8B cosh £ — sin (n + B) 
= — sinh ¢ tan~! ——_______—_———__ 
T cos B sinh & 


Let also sin 8 cosh §=p, cos B sinh E=q, n+8=y. Then 


1 ad — sin y 
Ri + 1S; f — sinh é tan“! L Eeetoved 5 e**ndy 
0 


T T q 


P sinh ff 7 ae sin y P sinh &e—**6 
———— tan- 
0 


I 


—_—_— e4dy = ——_—___ [,. 
q ™ 


9 
Tr 


Here the limits of integration need not be changed, since the integrand is a periodic 
function of period 27. To evaluate J;, use is made of integration by parts, with 





u = tan“! Le td dv = e**¥dy. 
Then : 
dy = — anne a , tae tooo ei*¥ (k ¥ 0) 
(p — sin y)? + ¢? 
and 


[2x Qe 
I, = uv} -f vdu. 
0 0 


But, since both u and » are periodic, their product evaluated over the period 27 is 


» } al cos wy d 
I, = -f vdu = — +f Leer oe eiky, 
0 kJo (p—siny)? + q’ 


zero; then 
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Replacing the sine and cosine by their exponential equivalents and transforming, one 


has 
2 oF (e?*¥ + 1)d(e¥ 
Sicai aa etk¥ —— jn da ) —_- 
kJo (ec? — 2ipe'¥ — 1)? — 4q7e?¥ 





This is easily seen to be a rational function of e*¥, the denominator of which, the 
difference of two squares, can be decomposed into two quadratic factors with rela- 
tively simple roots, so that the transformation by partial fractions can be used to 
obtain the following result: 


1 al 1 1 1 1 
: 2kdJ 9 e'¥ — efei® = pit — e-tgiB = gid 4 ete Bh gid 4 eke 8 (e"*) 


Thus the integral breaks down into four integrals of the form 


2n t* 
[gee 
0 t-—c 


where t=e'¥, and c isa complex constant of the form +e+*+*, Now if the indicated 
division of t* by t—c is performed, a quotient which is a polynomial in ¢ and a re- 
mainder c* result. The polynomial is integrated into another polynomial in t=e*, 
and the value of this second polynomial between the limits 0 and 27 is zero because 
of the periodicity of e‘¥. Thus the remaining terms are of the type 


Qn c*dt 2r ck 
<— =| ——de); 
0 t—c 0 ev —¢ 
namely, 


1 ml (ete) * (e~fe*) & (- efe— #8) * (- ee) * 


2k — eet 














I, = Jae, (A) 


ev — efet8 ev _— e fe ev + efe—# ev + ete ¥ 


The value of the resulting terms can be obtained more easily by considering the 
corresponding complex function of {=~+iw.* 


d(e*t) : etdt 
c* — = ict a aaa 
cet —¢ cet —¢ 


along the contour shown on Fig. 5. It is well known from the theory of the complex 
variable that the value of the above contour integral is zero if the pole of the integrand 
falls outside that contour, and is equal to 271Xic*X Res., where Res. is the residue 
of the integrand, if the pole lies inside the contour. 
Performing the integration around the contour, we obtain the following: 
(1) along the real axis —7 Sz, the complex integral reduces to the real 
integral to be evaluated (limits — and 7 are equivalent to 0 and 27); 
(2) the two integrals along the vertical paths,cancel each other; 
(3) the integral along w=A has a zero limit for A> ~. 





Therefore, 


* This treatment was indicated to the writer by his friend and colleague, Dr. H. G. Baerwald. 
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, O if the pole lies outside the strip 
2r c* 
f ——— d(e%) = {eZ 0 —rsysr 
ee vee iar ama lia ; 
| 2xi X ic* X Res. if it lies inside that strip. 
The pole occurs at e# —c=0, or {= —7 log c; ¢ is of the form 
c= ette® or c= — ette-# = ettei(rs), 
so that log c= +£+76 or log c= + §+i(7—8) and [=6 F72é or [=a7—B F 2. 
But the region x>0 corresponds to —<0, as can be seen from (2), and also 
0<8<7/2, so that the integrals whose pole has an imaginary part of the form +72, 


namely the second and fourth of (A), have the value zero. The poles of the first and 
third, on the contrary, fall inside the region of integration, so that their values are 


2mi X ic* X Res. 


It remains to evaluate the residue. This is found to be —7 by methods explained 
jn texts on the complex variable (Ref. 5). Thus the required integrals become 


2mei X ic* X (— i) = 2wic*. 


Thus, 
1 
I, = ok 2ri[(efe®)* — (— ete-#)*], 
and 
re... ert 
R, + iS, = —sinh §— [1 — e(r-#*], 
vs k 
Therefore 
P ekt 
R, = — (— 1)*—sinh £ sin kB —, 
T k 
P eké 
S, = — sinh  — [1 — (— 1)* cos hf]. 
T k 


When 6=0, the poles shown on Fig. 5 have real parts 0 and 7, respectively. In 
other words, one of the poles is on the contour itself. Besides, due to the periodicity 
of the integrand, a third pole appears with a real part equal to —7. This latter pole 
has, in general, a real part —a—{8, and is identical with the pole at r—. Thus, there 
are three poles in all, one wholly within the contour and two others, with equal 
residues, on the contour itself. Now it is easy to see that each of the latter contributes 
half its residue to the value of the integral, and since these residues are equal, the 
situation remains the same as if there were only two poles, both entirely within the 
contour, so that the case 8 =0 is not essentially different from 80, and it is sufficient 
to set @=0 in the above formulas for R; and S;. Thus 

(0) (0) aP (0) 


R, =0, S, =——sinhé for k odd, Si 
T 


= 0 for k even. 


Case k=0. For this case, the procedure is exactly the same up to the integration 
by parts. There, while u remains as before, dv=dy, v=y, so that J, becomes 
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= 29 tan“ p/q + qJo. 





— si |\2e Qn d: 
Io = y tan“ e-F +f qv — "7 
q 0 (p—siny)? + q? 
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To evaluate Jo, the same method as before is used, exponentials being introduced 
in place of the trigonometric functions: 





yas v(e?* + 1)d(e) 
Jo = 2i f , ; ° 
o (e — 2ipe® — 1)? — 4g%e?"¥ 


This can again be transformed into partial fractions: 


1 al 1 1 1 1 
Fe i —— ~ de). 
2 cos SJ. 5 —eie® ei — eteB eid + ete iB 7 ev + cal te") 


9 








Here we are dealing with integrals of the type 


f 2r yd(e'¥) if 2r pe'¥dy 
Saas tai pital : 
0 ev¥—e¢ 0 ev—e 


these can be treated, as before, by in- 
troducing the complex variable and 
integrating around the contour of ————— eee 4 
Fig. 5. Since the denominator of the 1 
integrand is the same as before, all 
that was said about the poles of the | 
partial integrals making up J; re- Y 
mains true. Therefore, the second and | 
fourth terms in the expression for Jo = 
contribute nothing, and the first and = 4 
third are each equal to 271 times the 

residue times constants. The residues, 

however, have here the value —log c, 














Fic. 5. 


so that 


i T 
ie oe ee a aa 
"pony i ri[—&— iB +E + i(e — B)] ae 28) 


and the imaginary term vanishes, as could be expected. Then 
I, = 2x tan— (tan 8 coth £) + w (# — 28) sinhé 


and 


P sinh & 
Ry = ——— [2 tan“ (tan 6 coth £) + (# — 28) sinh €]. 


T 
Half of this expression is the first term of the Fourier series: 
Ro FPsinhé 7 , 
— = ——— | tan“ (tan B coth &) + ry _ 8) sinh ‘|. 


2 T 
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For 8=0, there are again three poles, one on the imaginary axis and the two 
others with real parts +7. This location is, as in the general case k¥0, due to the 
periodicity of the denominator of the integrand. However, here the integrand as a 
whole is not periodic, so that the residues at the two poles on the contour are not 
equal, and the situation is not the same as for 8#0. The detailed computations show 
that Ro=0. 

Series for ¢://. Since the ratio 

p — sin (n + 6) 


o1/J:¢2/J = o1/¢2 
q 


is a simple trigonometric expression, the series for ¢1/J can be obtained from that 


for ¢2/J by term-by-term multiplication. 
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THE THERMAL-STRESS AND BODY-FORCE PROBLEMS OF 
THE INFINITE ORTHOTROPIC SOLID* 


BY 


G. F. CARRIER 
Cornell University 


1. Introduction. Elastic problems dealing with orthotropic materials have had 
considerable investigation in recent years, but up to the present time, such investiga- 
tion has been largely limited to a consideration of the problems involving thin plates 
of this material. 

In the present paper, two problems dealing with the stresses and displacements 
in an infinite elastic orthotropic solid are solved, and in each case the results are ob- 
tained in terms of three independent displacement potentials. The two solutions are: 
1) the displacement potentials arising from an arbitrary distribution of temperature 
within a finite region of the solid (the temperature being measured from an arbitrary 
datum) and 2) the potentials arising from an arbitrary distribution of body force 
within a finite region. Each of these problems reduces to the solution of three simul- 
taneous partial differential equations, which are transformed, through the use of 
Fourier integrals, into individual solutions for each potential. The expressions for 
these potentials are reduced to the form of Newtonian potential integrals for those 
cases where sufficient symmetry of the material properties exists to allow such a re- 
duction. In the more complicated cases, the results are still expressed in closed form 
in terms of definite integrals. 

2. The thermo-elastic problem. The conditions under which the thermo-elastic 
problem will be formulated and solved are the following. The material is to be homo- 
geneous, orthotropic, and elastic, throughout the infinite region, and is to be within 
that class of orthotropic materials which has three coefficients of temperature ex- 
pansion, a;, associated with the three principal directions of the material. The body 
forces will be taken as vanishing, since any problem involving both thermal and body 
force effects has a solution which is merely the superposition of the two individual 
solutions. The temperature distribution is to be an arbitrary function of position with 
the restrictions that this function must vanish everywhere outside some finite region, 
be continuous everywhere and be differentiable everywhere except on a finite number 
of surfaces. - 

The fundamental relations needed to formulate the problem mathematically are: 
the equations of equilibrium of an element of the material; the thermo-elastic equa- 
tions, that is, the relations between strains, stresses and temperature; and the rela- 
tions between strains and displacements. 

The equations of equilibrium are found by a consideration of the equilibrium of 
a rectangular parallelepiped of the material under general loading. Since these equa- — 
tions are independent of the type of material under consideration, they are given, as 


* Received Sept. 1, 1943. 
1 See, for example, A. E. Green, and G. I. Taylor, Stress distributions in aeolotropic plates, Proc. Roy. 


Soc. A 173, 163 (1939). 
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in the isotropic case for zero body force,? by three equations of the type, 


Oo - OT ay OT ss 
A+ —= = 0, (1) 
Ox oy Oz 








where the notation is the conventional one. 

The orthotropic material has been defined as one whose Hooke's law has the form 
indicated by equations (2), when T is identically zero. The effect of temperatures, 
different from datum, is to produce normal strains in the three principal directions 
of the material, as specified under the conditions of the problem. Hence, when the 
coordinate axes are taken parallel to the principal directions, the general formulas 
for the strains have the form, 


Cz = A102 + Ay20y + A130, + aT, +> | Yue = GUaaTyz * ** - (2) 
If we now define three displacement potentials, ¢;, such that 
_ 991 _ Ode ds 


’ v= ’ w= ’ 
ax dy dz 











u 


and such that @; and its derivatives vanish at infinity, the conventional definitions 
of the strains become, 


9 


Ou 0°64 Ov Ow 


e- = = is ¢e «= ¥ = _—_- = 
Ox dx? “as «Oy (Azdy 


($2 + $3), » Paths (3) 








Combining now, equations (1), (2), and (3), we obtain three equations of which 


the following is the first: 


0 0? 0? 0? 0? 0? OT 
=| (bu + ba nn a Sige Jon + on b+ on | = =~ * (4) 
dy? dy? dz? 


Ox Ox Vy 02? Ox 


Each of these may be integrated once to give,’ 





0? 0? 3? 0? 0? : 
fa + Oe; + Oe =)e + C12 — $2 + Cis—— 3 = —BiT,-++ (4a) 
Ox? dy? 02 dy? dz" 


- 


The arbitrary functions which appear in each of the foregoing integrations must 
each vanish, since, for example, in the first equation, all terms vanish when x is in- 
finite and y, 2 are finite, implying that all functions independent of x must vanish 
identically. 

Due to the convenient form of the boundary conditions, these equations are easily 
integrated by the following procedure. Multiply each equation through by e~ *(#t+ v"+4%) 
and integrate over the whole region, integrating by parts those terms containing 
derivatives of ¢;. This operation produces the following three equations, using the 
abbreviated forms defined below in equations (6). 


2 A. E. H. Love, A treatise on the mathematical theory of elasticity, Cambridge, 1934, p. 125. 
3 The b;;, cz; and 8; are combinations of elastic and thermal constants arising from the above opera- 
tion. The manner in which these constants appear in the second and third of these is easily deduced from 


equations (5). 
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(Diré? + Deon? + Desf?) Ei + Cion*E2 + Cisf*Es = BS, 
C12k? Ey + (beck? + boon?” + buf?) Ee + Co3E3 = B2S, (5) 
Cist?Ey + Cosn?E2 + (dsst? + Daan? + dah?) Es = B3S, 


E=- fff op je izttuntet dy dy dz, 


where, 


(6) 
= = fff Te ittettunt ab dy dy dz. 
Equations (5) are easily solved for the E;, and yield the expressions, 
E; = Fé, 9, OS, (7) 


where the F; become ratios of homogeneous polynomials in £, n? and ¢?. 
Noting now, that by their definitions, the E; are the Fourier transforms (in three 


dimensions) of the ¢;, we may write 


1 ce) 
eededacal” J J J F i(E, n, Set DS(E, m, )dE dn dk 
oT" — 


1 “ x 
_ rs fff F (En, Cetet ute dé dy df fff T(r, s, Lhe~*(ret rt D dr ds dt, 
oT —2 —o 


and the order of the indicated integrations may be changed to give, 


1 m 
¢; = 83 fff T(r, s, t)dr ds dt fff F,(é, 0, Cjeil(e—N E+ u-#) at (2—-OF dE dy df. (9) 
T° = —o : 


Since each F; (as defined by equation 7) is a ratio of second order polynomial in 
j y eq 
£*, n?, and ¢?, to one of third order, we may write, 
Rik; 
age came rag 
2p2pP2 
R?RiR? 


(8) 


F;=B 


where Ri =)zé*+ yin? +¢2, and where the Ax, ws, and B, are constants depending on 
the values of the constants appearing in the determinants defining F;, and hence, 
may be considered as known. Note that the j3, y?, for k=1, 2, 3, must be non-nega- 
tive, since no singularity may exist except at the origin. 

In many cases, the expressions for the F; may be reduced to the form, 


F;= >> ok (10) 
kai,2,3 Re - 

This will always be true when the problem involves a material which is isotropic in 
a certain plane (for example, a laminated plastic) unless identical values of R* recur 
in the denominator. This may be seen by noting that since the denominator of F; 
must be invariant under a rotation about the z axis due to this isotropy (the plane of 
isotropy is here taken as the x, y plane), &? and 7? must occur in the combination 
£+-?, and hence, A\x=mx, and the Rj become essentially binomials. The reduction 
of F; to the form of equation (10) is, in this case, merely a matter of evaluating A jx 
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When equation (10) does hold, the integration proceeds as follows: Using the con- 
ventional vector notation and the new coordinates with the subscript k (where 
Ee=Axnk, x, =xXAT!, re =Phe'!, etc. and where Mm, =i(x,—7rx) +)(ye—Sz) +h(2e—Ex)), the 
integral over &, y, and ¢, of equation (9) defining Green’s function G, may be written, 


; .~ d& d 
G,(x, y, 2,7, 5,8) = Sfp. ‘. — gimk Rk brid nd (11) 
k= a R At Mk 


If we now change to a spherical coordinate system in which y is the angle between 
m, and R, and 6 is the polar angle about m,, this integral becomes, 


w= fff x —  gimeRe cos y sin y dy dé dR,, 
om 


where the integration now takes place over, OS y Sm, 0S 5S 27, OST Ri < ~. The ele- 
mentary integrations over y and 6 produce 
; An. lik °@ sin mR; 
G;= b> — —— dR,, 
k=1,2,3 ArmmE Ry 
which is known to have the value, 
. 2 Aj, 1 
G; = 27? — —-- 


k=1,2,3 Nek Mi. 


Now transforming the remaining terms of equation (9) to the coordinates with the 
subscript k, and substituting the above value for Green’s function, we obtain, 


” i », ogee Eth rn Shc a 


Hence, the problem, wherein T(x, y, z) represents the temperature distribution, 
becomes the problem of evaluating the Newtonian potential function corresponding 





to a mass distribution of, 


Ask 
p = —T(AcXk, MeVey Zk). 
4r 


For an isotropic material, the ¢, become alike, and are given by,‘ 


a ito Tr, s, t) 
¢, = tp fff — a dr ds dt. 
4n 1 —9 We rier rs 


In the evaluation of Green’s function for those cases where the denominator of F; 
has a multiple root, it is convenient to introduce the notation 








g? 0? 0? 


= > Gin, By Ome e's} ° 
: ax? | dy? as? 





a 


In this case, » integrals of the form, 


‘J. N. Goodier, On the integration of the thermo-elastic equations, Phil. Mag. (7), 23, 1017 (1937). 
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= AiR - 
Gin _ fff a et l(a E+ (ys) at (2-H 1 gE dn dt (13) 


2 


must be evaluated, provided Rj R}. Using the above notation, the equivalence of 
the following to equation (13) may be verified by substitution: 
A in 
AGin = af ff aie etl(a—n) E+ (yma) at tO dE dn dt. (13a) 
R; 
The integral involved in this equation is, however, the same as that appearing in 
equation (11), so (13a) becomes, 
Cn 
AG in = At i Se a — ’ 
V (x2 — 72)? + (v2 — $2)? + (2 — te)? 


where Cj, is an easily evaluated constant. Substitution will again show that, 








Gin = C jnAv/(%2 — 172)? + (ye — 52)? + (22 — é2)’, 


is equivalent to the above equation, and hence the ¢; are given by, 


1 . 
o;= 8x3 fff T(r, 5, )Gin(7, 5, t, x, y, 2)dr ds dt. (14) 
TT on 


In those cases where F; cannot be reduced to one of the foregoing convenient 
forms, G; is more difficult to evaluate. Since no explicit form has been found for this 
function, other than complicated definite integrals, it is believed best to leave it in 
the form defined by equation (9). 

3. The body-force problem. As in dealing with isotropic materials, the solution 
of the body force problem may be shown to reduce to a form analogous to that of the 
thermo-elastic problem. To show this, we shall consider only the problem where the 
body force is directed parallel to the x axis, noting that the general solution is ob- 
tained by the superposition of three such problems. 

Equations (1) and (2) are modified to contain the body-force function, X, and to 
eliminate the temperature terms. Equations (4) are then obtained again, where now 
the right hand sides are replaced respectively by, X, 0, and 0. 

The ¢; will not, in general, vanish at infinity in this problem, hence the procedure 
needs a slight modification. The second and third of these equations are integrated 
with respect to y and z respectively and then differentiated with respect to x. This 
yields equations (4a), where again, X, 0, 0, appear on the right and where the @¢; are 
replaced by 0¢;/8X. The procedure is now identical with that of the thermal problem, 
and the ¢; are found by the expressions analogous to equation (14). 

4. The two-dimensional problem. If we carry through in two dimensions the 
procedure used in the previous sections of this paper, we arrive at an equation which 
is identical to equation (8) except that z, t, and {, no longer appear. The expressions 
for F; are now simpler in form, being given by, 


7 Nie? + win? 
"(NEP + win?) AE? + win") 
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which may always be reduced to the form 


F;= > 


k=1,2 R? 


A jk 





unless Re=R3. (R2=)z +n’). 
Before changing the order of integration, we differentiate equation (8) with re- 
spect to y. The integral form of Green’s function becomes then 


OG jx 3 ine* (=e un) 
——_ ff A ei ages gd dé dn, (15a) 
Oy —0 ze? + pin? 
unless R2=R3, in which case, 
aG; ines (zttvn) 
A, — = a ff A; ——— dé dn. (15b) 
oy AoE? + pon? 


This latter expression is, of course, derived by the same reasoning used in the three 


dimensional problem. 
Equation (15a), after the introduction of the coordinates with the subscript k, 


can be written in the iterated integral form, 


£% 
OG 5% oe fA ik, os COS NkXk Nk & 
- = —— sin i Yu One - ———— 6 { — 
oy 0 AE 0 £, Nk 
1+ 


which is known to be equivalent to, 


OG jx - 4A jk TT : 
——_—_—_- = — sin Venk* e~!knkl dn, 
oy 0 AK 2 











and this integral yields, 
OG jx 20A jx Vk 


dy + 9 


7 








or 
WA jk 
+ 2 2 
Gy, = —— In (xk + yi), (16) 
Agu 
and we obtain the familiar two-dimensional logarithmic potential. 
Equation (15b), then becomes, in an analogous manner, 





TA i 2 9 
AxG; = —— A, In (x2 + y2) 
AkMe 
or, 
WA ; “ re 2 ; 
G;= Ao[ (x2 + y3) In (x3 + 43) ]. (17) 
Q2rkME 


Hence, Green’s functions are determined for each two-dimensional problem involving 
thermal stress or body forces in the infinite plate. The usual methods of superimposing 
plane stress (or strain) solutions may be utilized, of course, to solve the corresponding 


problems for the finite body. 


STRESSES IN THE DIAPHRAGMS OF DIAPHRAGM-PUMPS* 


BY 


A. M. BINNIE 
Engineering Laboratory, Oxford 


1. Introduction. Complete prevention of leakage from a reciprocating pump is 
difficult to ensure over a long period of working. When the fluid to be pumped is 
of such a nature that no leakage whatever is permissible, some modification in the 
design of the pump is essential, and under these circumstances a diaphragm-pump 
may conveniently be used. This in its essentials consists of two chambers (Fig. 1) 
attached to a modified reciprocating pump. The chambers are of conical form rounded 
off at the apex and at the base, and between them a diaphragm is clamped at its edge. 
For high-pressure operation the diaphragm is a very thin steel disc. The fluid to be 
pumped passes through one chamber, connexion to inlet and exhaust valves being 
made by means of a number of small ports. The other chamber is connected similarly 
to a single-acting reciprocating pump, which is not fitted with valves. This chamber, 
the pump cylinder and their connecting ports are filled with a liquid (commonly oil), 
and thus motion of the piston of the reciprocating pump causes the diaphragm to be 
pressed alternately against both conical surfaces, thereby producing the desired 
pumping action. The inevitable leakage of oil past the piston is made good by means 
of an auxiliary pump. 

An approximate method of calculating the stresses in the diaphragm is explained 
below, hence the size of the chambers may be so designed that the fatigue strength 
of the diaphragm is not exceeded. In section 2 the deflexion of the diaphragm is 
taken as sinusoidal, in section 3 as a cubic, and in section 4 as following a Bessel- 
function relation Attention is confined to the stresses which result from distortion 
into the same shape as the chamber, no regard being paid to the local stresses round 
the ports. 

2. Stresses when the transverse displacement is sinusoidal. In general the dis- 
placement of the diaphragm from its unstrained position has not only a transverse 
but also a radial component; therefore it 
does not seem possible (except by relax- atti 
ation methods) to calculate the stresses (——_i_,_if = 
for a specified shape of chamber. It is 
necessary to assume a reasonable ex- 
pression for the transverse displacement Fyg. 1, Arrangement of diaphragm and chambers. 
w, from which the corresponding radial 
displacement u and the stresses will be obtained; and, when both wand w are known, 
the shape of the chamber is determined. 

With the axes shown in Fig. 1 we shall in this section take w as specified by 


Wo wr . 
w= “(1+ 00s), (1) 
2 a 
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where a is the radius of the diaphragm and wo the maximum value of w. This expres- 
sion satisfies the conditions that the slope must vanish at the centre and at the edge. 


For u we assume as an approximation! that 
u = r(a — r)(Cy + Cor), (2) 


where C; and C; are constants which will be determined by the principle of minimum 
strain energy. The conditions that u is zero at the centre and at the edge are auto- 
matically fulfilled. Now the transverse displacements are many times the thickness 
of the diaphragm, hence large-deflexion theory must be employed. In Timoshenko’s 
notation (loc. cit.) the radial and tangential strains are thus 


me (=) 

eo —4—i-—}, 

dr 2 \dr 
u 


= * 
r 


(3) 


The diaphragm being very thin in comparison with its radius, the strain energy in 
it due to bending may be neglected in comparison with that due to the stretching of 
its middle plane. Hence the strain energy in the diaphragm is 

wEh 


Vi= f r(ez + e2 + 2ve,e,)dr. (4) 
1 - yp? 0 





Here E denotes Young's modulus, v Poisson’s ratio, and & the uniform thickness of 
the diaphragm. On putting (1) and (2) into (3) and inserting the results in (4) we 


obtain 














wEh [1 . 3 7 an 
vy,;=>— l— Cjat + —C,C.a® + — Coa’ + Wo 
1—v*L4 10 60 1024 
wa : 7 a 1? 1 5 a | 1 
+fa(i-=+(+—)) +cm(=-5+(4+-)) . (5) 
8 6 6 2 + iz 12 + 
Now 
OV, dV, 
re] Oi OC. 
hence 
C 25 we {= 1 17 (= 4 1) 
me's NS ; a 
15wi (x? 1 
C= - {= +13 - »(= +1) 
128a* | 3 3 ‘ 
In the remaining calculations we will consider the case vy=0.3 when (7) reduces to 
ws wo 
C, = 1.06—,» C;: = —1.76—- (8) 
a’ a 


1 Cf. S. Timoshenko, Theory of Plates and Shells, McGraw-Hill, New York and London, 1940, chap. 
IX. 
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The radial and tangential tensile stresses are 








E 
(e+ ve), of= (e, + ve,). (9) 


o; = 
1 — »? 1 — »? 


Then, with the aid of (3) and (8), (9) becomes 





Euj r r _. - 
o; = {1.st— 7.11% + 6.385 + 1.36 sint 7, 
a 





a? a a (10) 
Ew} r r — | 

o; = *)1.51 — 4.95 —- + 3.67 — + 0.41 sin? at 
a? a a? a 


These expressions are plotted in Fig. 2, from which it will be seen that the maximum 
stress occurs at the centre and is given by 1.51 Eug/a*. If wo/a=1/35 and E=13000 
tons/sq. in., the maximum stress is 17 tons/sq. in., which for a good quality steel is 
a reasonable working stress. Finally we will examine the shape of the chamber corre- 
sponding to (10). It will be noticed from (2) that u is zero not only at the centre and 
at the edge but also at r/a= —C,/(C.a) = 0.60. On differentiating (2) it appears that u 
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Fic. 2. Radial and tangential tensile stresses in the diaphragm. 
— —-—— sinusoidal displacement. 
—:‘—-:— _ cubic displacement. 
—:+-—-+-+— Bessel-function displacement. 


has a maximum value 0.12u3?/a at r/a=0.24 and a minimum value —0.06w6/a at 
r/a=0.82. Hence the greatest radial difference between the shape of the chamber 
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and w as given by (1) is 0.12 w%/a, which for wo/a =1/35 is only 0.98 X10—a. For a 
chamber of normal diameter this difference is too small to be appreciable in manu- 


facture. 
It is of interest to determine the strain undergone by a radius of the diaphragm, 


and for this purpose an accurate method of rectification is available. For a chamber 
of sinusoidal shape a radius is extended to a length s given by 


2a re ° mr\ 112 wl 
= ) paar f (1 — p* sin? ) i(~), (11) 
~ 0 , a a 


: rw A rws 
where jkr re i+ ; 


Since p is small, the first bracket in (11) may be expanded by the binomial theorem 
and the complete elliptic integral replaced by 


s-a p* 139* 
me) ay Tse a Se (12) 
a + 64 


For wo/a = 1/35 this strain amounts to 0.05%, hence in a steel wire distorted into this 
sinusoidal form the tensile stress would be only 0.0005 K 13000 =6.5 tons/sq.in. 

3. Stresses when the transverse displacement follows a cubic relation. To esti- 
mate how far the stresses depend on the expression assumed for w, we will in this 
section replace (1) by the cubic 


3r? 2r8 
w= w(1-— +=), (13) 


a? a* 


This equation satisfies the same four boundary conditions as (1), and the greatest 
difference between the two is approximately 0.010w» at r/a=0.28 and 0.72. After 
employing (2), (3), (4) and (6) we find that 





wEh 1 3 a. 9 wo 
ae E Cia‘ + —C Cea + — Cia® -— — 
i-0Ls 10 60 35 a? 
3 
+ = wha{2Ci(3» — 1) + Cza(3v + 1)] : (14) 
é at 
c,= 3% 23-15), 98 41 — 5) (15) 
= i ov), ,=> =— = Die 
' 36a3 ‘ 56a! 
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For v=0.3, (15) reduces to 
2 2 


Wo Wo 
C; = 0.99 —» C; = = 1.62—, (16) 
a’ a‘ 


and the stresses obtained from (3) and (9) are 


ll 


Ew? r r? r r! 
Or — {1.42 — 6.61 — + 25.67 — — 39.56 — + 19.78 —>, 
a? a a? a’ a‘ (17) 
Ew} r r? r 
o1=- 1.42 — 4.60 — + 9.32 — — 11.87 — + 5.93 —>. 
2 a a? a’ a‘ 


a* 





From Fig. 2, in which these expressions also are shown, it will be seen that the maxi- 
mum stress is slightly smaller than that obtained in section 2. 

4. Stresses when the transverse displacement follows a Bessel-function relation. 
Lastly we will take w as given by 


w= Wol{Jo(kr) — m}, (18) 
where k=a/a, a=3.83 +--+ being the first positive root of Ji(x) =0, 
m = Ja) = — 0.402---, 


and 


Wo = wo/{Jo(0) — Jo(a)} = wo/1.402---. 


This equation satisfies the four boundary conditions, and it gives a displacement 
which, unlike those previously considered, is unsymmetrical about the line w=w/2. 
Except at r=0 and r=a the displacement is everywhere less than that specified by 
(1), the greatest difference between the two being approximately 0.019 wy at r/a = 0.53. 
The same procedure as before leads to 











wEh [1 3 7 Wok ¢* 
Vyi,=— | Cia* + — C,C.a® + — Cia® + f rJi(kr)dr 
1— v*L4 10 60 0 
ky 372 ka 
+ wif (Cra + 3C)) 
12 
3 a 
+ r (2 + v)(C2a — C) (f Ji(kr)dr — asa) |, (19) 
5Wo : 
C,= —- | 6¢2 + m4 f Ji( kr)dr — aJ%(ha) — pra Fi(ba) |, 
a 
sw? ~ (20) 
C.= - , — | 9 + rf Ji(kr)dr — aJi(ha) o brea) |. 
a 0 


If we take v=0.3 and f[*Jo(x)dx =1.2599,3 (20) reduces to 





* After some fruitless attempts to evaluate this integral, I asked Professor G. N. Watson whether 
it was expressible in any simple form; his reply was that he thought not, and he computed its value to 
15 places of decimals, his result being 1.25990 97359 05768. The value 1.2599 is sufficiently accurate for 
our present purpose, 








42 A. M. BINNIE 
we wi 
C, = 1.01—, C, = —1.74— (21) 
a’ a‘ 


and the stresses are 


I 


Ew? oe, r? r 
oy = —-21.44 — 6.95 — + 6.31 — + 4.10/2(kr) $, 
a a a? (22) 
Eu} r r’ " 
o, = ——41.44 — 4.84 — + 3.63 — + 1.23Ji(kr) 2. 
a? a a? 
These stress distributions, which are plotted in Fig. 2, are in close accord with the 
results obtained in sections 2 and 3. 
5. Conclusions. The following conclusions emerge from the above calculations :— 
(t) For the three kinds of displacement considered, the maximum stress in the dia- 
phragm is at the centre and is about 1.5 Ewé/a?. 
(tt) The stress distributions due to the three kinds of displacement do not differ 
widely. Hence, if it is decided to use one kind, and small errors are made in the 
difficult process of machining the chambers, no great alteration in the stresses 


will result. 
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THE INTRINSIC THEORY OF THIN SHELLS AND PLATES 
PART II.—APPLICATION TO THIN PLATES* 


BY 
WEI-ZANG CHIEN 
Department of Applied Mathematics, University of Toronto 


7. The general equations for a thin plate. We shall now investigate the equations 
of equilibrium and compatibility for a thin plate, not necessarily of constant thick- 
ness. First, we shall introduce the condition that the system is a plate, i.e., its middle 
surface in the unstrained state is plane. We have therefore 


bes = 0, Resy = 0. (7.1) 


Furthermore, in order to simplify the problem, we assume in the following sections 
that the body force forms a parallel vector field, and therefore (3.38), (3.39) are 
satisfied ; this is true for most practical problems. 

Substituting (7.1) and the conditions on body force into (6.34) and (6.35), we have 
three equations of equilibrium for a thin plate 


= 2A @oyPah + AGS (Qerh®) joy + ART rwGrsQoyh® i 
a 


o 
+ 2Xjo)h + (Oth) |. + yng a ®qnQ°h = Olea)» (7. 2a) 


2A (Perh) |p — AGS” (Grugnsh*) |p + $a" QeyAt) (Qrsh*)), + Po 
a ei a 


o 
+ 2Xigh + a a*?(0°h) |, + (a™¥a2? + 2a%a™)0,qr3h = Of), (7.2b) 
— o a 


where the O-symbols have the following magnitudes, 


Oks); X°ph, X°Oh, p’gh, O*gh, Oagph, Pqh?, X°h’, On}, qgPh', gph’, qX h', gh’. (7. 3a) 


Oty; pth, Oph, Xph, XOnh, 02h, Pqh?, On’, X hi’, ph’, gh. (7.3b) 
We recall that 
1 
Aig? = {oatta® + (1 — a)astah}, (7.3¢) 
nail 
2(2¢ — 1) 
Agen — ——_* gin gabry 4. ahs goby, (7.3d) 
3(1 — oa) 


Similarly, substituting (7.1) into (6.43) and (6.44), we have three equations of com- 
patibility for a thin plate 
2niiQasiy(1 + 2nijin\PaPx + 2a™p-a) 
a 
a 2nioigsx(a”™ + 2niinio}Pps) (Parl + Prie — Pay\n) = 0, (7.4a) 
a a a 


* Received June 12, 1943. Part I of this paper appeared in this Quarterly, 1, 297-327 (1944). 
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B By mr 
(2nfpinio|Por 1 t nfo} VovQaa) (1 + 2njpNo}PrP ro + 2a**p,:) 
a 


pa By/ TA. w | Tw ~ 
+ nfoittio} (2njojnio}Prs + a®*)(Dspiy + Pyrip — Prp\r)(Pawis + Piola — PoBiw) = 0. (7.4b) 
a a a a 


H a 
The macroscopic tensors in (6.29), (6.30), (6.31) can be written as 


Co 





Te? = 2AG Pah — AG) aeqeh?+ —— av0°h + Oy, (7. 5a) 
—¢ 
Le = 3nijar,Agy ash? + Of), (7. 5b) 


T° = 2A (qysh3)|, + Oh + (a P* + a**P)q,,h? + 0%, (7.5c) 


where 
O78, = O78)(p2h, 02h, Oph, phi, Qh’, Xh', g*h'), (7. 6a) 
(7. 6b) 


Of) = Ot (aph’, X°h’, Qh, g*h’, Xqh'), 
O72, = O78,(Oph, 02h, OX h?, Ogh?, Pph?, PQh?, gph’, X°h*, Qh’, Xqh', g*h’). (7. 6c) 
Equations (7.2a, b) and (7.4a, b) are the six differential equations of a thin plate 
in the six unknowns Pag and qas. The next step is to introduce certain systematic 
approximations based upon the thinness of the plate, so as to obtain a set of differ- 
ential equations in simpler form. 

8. Classification of all thin plate problems. We consider a family of ~! thin plates 
of the same material, having an identical middle surface S¢ in the unstrained state, 
but different thicknesses; each is subject to the action of (i) external force systems 
applied at the edges, (ii) surface loadings on its two boundary surfaces, and (iii) uni- 
form body force throughout the plate. (This includes gravity, but excludes a centrif- 
ugal field.) We attach to the middle surface of each plate the same system of 
coordinates x*%, so that the fundamental tensor aag is the same for all plates in this 
family. We assign to each plate a value of a parameter e€, so that the thickness of all 


the plates can be represented by 
2h = 2ch(x', x’), (8.1) 


where 0<e<e, and the function h is the same for all the plates; for thin plates, e: is 
supposed to be small, but the basic idea of the method is that we seek solutions valid 
for all € in the range 0<e€<&. 

Equation (8.1) implies that the derivatives of the thickness at any point are of 
the same order of magnitude as the thickness itself. We shall call these plates “regular 
plates.” On the other hand, if the thickness and its derivatives are of different orders 


of magnitude, we have an “irregular plate.” The following theory is limited to regular 


plates only. 

We may suppose € chosen equal to the ratio of the average thickness to a selected 
lateral dimension (usually the smallest lateral dimension) of the plate. For a circular 
plate, € is the ratio of the average thickness to the diameter of the plate. For a rec- 
tangular plate, it may be chosen equal to the ratio of the average thickness to the 
length of the shorter side. 

It is important to observe that € is the only parameter involved. Except the funda- 
mental tensor agg and Poisson’s ratio ¢, all the other quantities occurring are func- 
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tions of €, and no quantity is “small” unless it tends to zero with e. (Young’s modulus 
does not appear, on account of the use of reduced stresses and body forces.) Thus for 
any “small quantity” V, we must have 


li = 0. ; 
lim W 0 (8.2) 


In order that a problem may belong to the theory of small strain, e;; must be a small 
quantity, and therefore 


lim e;; = 0. (8.3) 
e—0 


It follows that pag must also be a “small quantity,” depending on ¢ like V in (8.2). 
But this is not necessarily true for qas. 

It is understood that all conditions (such as reduced edge forces, reduced surface 
loadings, and reduced body forces) depend on €¢ in such a way that (8.3) holds. We 
shall assume that Q‘, P‘, Xj9 vanish at least as fast as €, and are in fact power series 
in ¢. This assumption implies that the derivatives of any of these quantities with 
respect to x* are of the same order of magnitude as (or higher order of magnitude 
than) the quantity itself. Hence we write 


Q° ai +R Q%se", Qe ba > Ove’, (8. 4a) 
s=kg s=k 

Pe ve pe Pre’, Pe = Zz P%,,€*, (8. 4b) 

Xo) = Zo Xe t01€", Xo) = } X%s) (01€* (8. 4c) 
8=)9 s=] 


oe . i . 
where k, ko, , 0, j, jo are integers greater than zero, and P%,), Qf), X{.yo) are functions 


of x*, independent of e. 
Similarly we assume that the traction, shearing force and bending moment ap- 
plied on the edge curve can be represented by 


s=l 


Te = > The, (8.5a) Le = > Le, (8.5b) To = > Toe, (8.5c) 
s=t s=u 

where ¢t, u, / are positive integers, and rr. ES. ay are functions of position on the 

edge curve, independent of e. 

Now the problem is to find the behaviour of the family of ©! thin plates under 
the action of a given family of external force systems (8.4), (8.5). Given an external 
force system defined by (8.4), (8.5), we seek solutions of the equations of equilibrium 
(7.2) and the equations of compatibility (7.4) of the form 


Pos = >, Pooast’, (8. 6a) Gas = >, Qarase’ (8. 6b) 


&=p s=q 


where p and g are zero or positive integers, and P(.)ag and q (sas are functions of x‘, 
independent of ¢. Only those problems admitting solutions with p>0 belong to the 
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theory of small strain. On the other hand, g may be zero; then we are dealing with a 
finite deflection problem. 

The usual discussion of plate theory is based on the deflection, i.e. the normal 
displacement of a particle on the middle surface. The present method is intrinsic, and 
the general equations contain no explicit reference to the displacement. However, 
since @as corresponds to change of curvature (i.e. curvature of the middle surface 
after strain), it is clear that finite values of qag correspond to finite deflection and small 
values of gas to small deflection. Similar remarks apply in the case of shells. Hence, 
in classification, we may use the familiar word “deflection” when referring to the order 
of magnitude of qag. 

The assumed forms (8.6a, b) imply that the derivatives of pas, qas with respect 
to x* are of the same order of magnitude as the quantities themselves, or of higher 
order. In fact, pag and qag expressed by (8.6a, b) represent the behaviour of the 
family of ©! thin plates under the action of the given family of P‘, Q', Xj, 
T«8, T«0, [8 defined by the equations (8.4), (8.5). It is understood that if P#, 03, X op 
T«8, To, [8 are identically equal to zero (i.e., k, ko, m, No, J, jo, t, uy L= oo), then 
Pag and gag vanish (i.e., p, g= © ) everywhere; this corresponds to the unstrained state 
of the plate. This means that self-strained plates are not discussed. 

In a thin plate problem, we are to regard the numbers k, ko, m, mo, j, jo, t, u, 1 as 
given; the initial step towards solution would appear to be the determination of p 
and q, for then we could simplify the equations of equilibrium and compatibility in 
the first approximation by picking out the principal terms in ¢ from equations 
(7.2a, b), (7.4a, b). But owing to the partial indeterminacy of p and g, this method 
is not successful. 

It is much simpler to solve the problem in the reverse order. First we assign in- 
tegral values to p and q. The values of k, Ro, , mo, 7, jo are fixed by the conditions that 
Xewtoy XCyo}r Pony» Pon)» Oey» OG) should contribute to the principal parts of 
(7.2a, b), without dominating these equations to the exclusion of pag and qas. The 
equations of equilibrium and compatibility in the first approximation are then ob- 
tained by picking out the principal terms in ¢ from equations (7.2a, b), (7.4a, b). Then 
the values of ¢, u, ] are automatically fixed through the expressions (7.5). 

We shall now discuss the classification of thin plate problems based on assigned 
values of p and q, so that the principal parts of (7.2a, b), (7.4a, b) in the first approxi- 
mation are different for different “Types.” The classification is shown graphically in 
Fig. 3, where permissible pairs of (p, g)-values are represented by circles. As indicated 
in (8.6a, b), we consider only non-negative integral values of p and g. Since, however, 
p=0 corresponds to finite extension of the middle surface, we must omit the (p, q)- 
points on the g-axis. 

It is found that the points in the (p, g)-plane break up into twelve groups depend- 
ing on their positions relative to the division lines AD, AB, OC and the p-axis. For 
any point (except g=0) on the line AD, it is easily seen from inspection of (7.2a) 
that the first and second terms are of the same order of magnitude and prevail over 
all the other terms, with possible exception of those involving Xjq, P*, Q*. For any 
(p, g)-point (except g=0) above AD, the second term in (7.2a) dominates, and for 
any (p, g)-point below AD, the first term dominates. For the point A, the first three 
terms in (7.2a) are of the same order of magnitude and prevail over the right hand 
side. For any point on the p-axis above A, the second and third terms in (7.2a) are 
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of the same order of magnitude and prevail over the other terms. Thus the principal 
part of (7.2a) takes five different forms depending on the position of the (p, q)-point 
relative to the line AD and the p-axis. 

Similarly, the form of the dominant part of (7.2b) depends on the position of the 
(p, g)-point relative to the line AB and the p-axis. Finally, the form of the dominant 
part of (7.4b) depends on the position of the (p, g)-point relative to the line OC and 
the p-axis. The equation (7.4a) has no division line, since the term nfo} VaB | 7 dominates 
for any position of the (p, q)-point. 
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Fic. 3. Classification of thin plate problems. 
p=order of extension of middle surface, g=order of change of curvature of middle surface. 
(Type P12 is not indicated in the diagram, since for these problems, g= ©, and consequently the corre- 
sponding points lie at infinity to the right hand side.) 


It follows that the (p, ¢g)-plane is divided into twelve regions, so far as permissible 
non-negative integral values of p and g are concerned, and so the complete classifica- 
tion of all thin plate problems involves consideration of twelve types (Types P1—P12). 
Type P12 is not indicated in the diagram, since for these problems, g= #, and con- 
sequently the corresponding points lie at infinity to the right hand side. 

Although the classification gives twelve types, four of these (Types P3, P6, P7, 
P8) are less important than the others. They represent overdetermined problems, in 
which the number of equations exceed the number of unknowns. Such cases can occur 
only when very special relations connect the body forces and surface forces. 

These twelve types may be described as follows: 

(1) Problems of finite deflection (¢=0), Types P1—P3. 
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(2) Problems of small deflection (q21, p=1; gq=1, p=2; g21, p>2q), Types 
P4-P8. 

(3) Problems of very small deflection (¢22, 2g=p22), Types P9-P11. 

(4) Problems of zero deflection (¢= «), Type P12. 

In order to save space, we shall not discuss all the twelve types in detail. The 
discussion of Types P1, P2, P3 will serve as an example. The results for all types 
are summarized in the tables in the Appendices at the end of this paper. The prin- 
cipal parts of the equations of equilibrium and compatibility are shown in Table I, 
and the orders of magnitude of the external forces and the principal parts of the 
macroscopic tensors in Table II. 

It should be noted that the theory of generalized plane stress [1, 2], the Lagrange- 
Kirchhoff theory of small deflection [3, 4, 5], and the von Karman theory of “large” 
deflection [6] can be derived respectively from the Types P12, P11, PS. 

We shall devote the next section to discussing the problems of finite deflection 
(P1—P3). All results for these types are new, and may prove particularly interesting. 

9. Problems of finite deflection (¢=0), types P1— P3. 


(a) Type P1: g=0, p=1. FINITE DEFLECTION WITH DOMINANT EXTENSION IN 
THE MIDDLE SURFACE 
General equations. By the condition that, in the first approximation, (7.2a, b) re- 
ceive significant contributions from P?,,), Pt), XG.joyr X Gyo» Oley» OG» WE must have 
No = n = 2, ign =j = 1, k=k=1. (9.1) 
Therefore, we obtain from (6.23) 


hi) = he + O(e*), hy_) = he + O(€); (9.2) 


consequently, the common assumption that the middle surface of the unstrained 
plate is deformed into the middle surface of the strained state is justified in the first 


approximation. 

We now substitute (8.1), (8.4)—(8.6) into (7.2), (7.4). The lowest power of € oc- 
curring is €? in (7.2), and e€® in (7.4). The corresponding coefficients give rise to equa- 
tions of equilibrium and compatibility in the first approximation as follows: 


J E z 1 — 2 ve 
— 2A a c0yerPayerk + Ply + 2X Gyo + Qtyh) ie + — a? qonQiyk = 0, (9.3a) 
a — 


~ rs o m 
2AM (p21) 2rF) |p + Poy + 2XGiy4 + 1 wee a*(O8h) |p 


+ (a®a°7 + 2a®%a™)gioynQiyrhk = 0,(9.3b) 


nj (oar = 0, (9.3c) 


nig} Mo} (0) 08 (0a = 0. (9. 3d) 


We may remark that all quantities in the above equations are finite, i.e. independent 
of e. The macroscopic tensors in (7.5) can be written as 
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a a 
T°! = (2AGy*Payah + tus a%®Qi,h)e + O(e), (9. 4a) 


L*? = 3niawAty Gonski + O(¢), (9. 4b) 
T2 = Of yhe? + O(e*) if Qh) x 0, 
T? = {Oh + ZA (Qoonsh®)|}@ +O) if Q4, = 0. (9. 4c) 


Equations (9.3a, b, c, d) form a set of six equations for the six unknowns q(o)ag and 
Pajas- From (9.3d), we see that the total curvature of the middle surface does not 
change for the first approximation. Consequently the strained middle surface in this 
type of problem is to be regarded as a developable surface. 

Equations (9.3a, b, c, d) can be further simplified. Since for a plate, R,as,=0, the 
order of the operations of covariant differentiation is immaterial; consequently, from 
(9.3c), we have 

J (0a = Wo) | a8- (9.5) 


Here w@) is an unknown function of x*, which satisfies, in consequence of (9.3d), 


njoiMioj (0) |pdW (0) ay = 0. (9.6) 
a a 


The existence of wo), satisfying (9.5), is easily proved by temporary use of special 
coordinates (rectangular Cartesians). The last equation is, in fact, the famous differ- 
ential equation [7] of a developable surface in the curvilinear coordinate system, and 
Wo) may be called the deflection function. If (9.6) is satisfied, qi) xy is given by (9.5). 
There still remain the three equations (9.3a, b) for the three unknowns puy,sy. We 
can handle the problem indirectly by means of T¢}. This is the coefficient of the lowest 
power, €?, in the series for 7**, and by (9.4a) 





- o » 
To = 2A Pp 1) aah + Fa a%®Oiyh. (9. 7) 


1— 


We note that this is a symmetrical tensor, so that it has only three independent com- 
ponents. Substituting (9.5), (9.7) into (9.3a, b), we have 


_ Te) W (0) |r + PO) }- 2X (iyo) + (QiyA) | + a wo) | Qik = 0, (9. 8a) 
Te)\2 + Pe) fe 2X%iyo) # + (a™ae? + 2a™a*?)w (oj 2Qc) yh = (), (9.8b) 


To sum up, for problems of type P1, we have a set of four equations (9.6), (9.8a, b) 
in the four unknowns, wo) and 7%. 
Special case. The following special case is interesting. If 


Po) = Xi = Oy = 9, (9.9) 
then by (9.8b) there exists an Airy function x, so that 


TE) = Mop MloiX 21 =- (9.10) 
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This is easily proved by temporary use of special coordinates (rectangular Cartesians 
[2]). Consequently, (9.8a) can be reduced to the form 


x nd ji SS £ 
— momo} X (2) 1a (0) 10x + Play + 2XCiyohk + a? woop nQtyk = 0. (9.11) 
a a a 


The problem is now to find x2) and wo) as functions of x* satisfying the two non- 
linear partial differential equations (9.6) and (9.11). In rectangular Cartesian co- 
ordinates, the equations (9.6) and (9.11) may be written as 


W (0),12W(0),12 — W0),11W(0),22 = 9, (9. 12a) 


o 7 0 -0 - 
2x (2),12W (0),12 — X(2),11(0),22 — X(2),22W(0),11 + Oey hAw oy + Pox» +2X (1) [0)4 = 0, (9.12b) 


where the comma indicates partial differentiation with respect to x*, and A is the two- 
dimensional Laplace operator. The macroscopic tensors are given by 





T" = x 2),22e” + O(e'), T?? = x (2),11€" + Oe), ! 
; oa (9. 13a) 
T!? = JT) = — x 2),12€7 + Of), 
pu = [** = D1 ait o)W (0),12€° + NM) 
L® = D(woo),11 + oW.o),22)€? + O(e), (9.13b) 
L* = — D(wo),22 + oWo,1)e + O(e), 
7* = ‘(1 — ¢)(Dwo),12),2 + (Dw), + oDw 0),22),1}€ + O(e),) (9. 13c) 
> eee 
T?? = ‘(1 - o) (Dw (o),12).1 + (Dw 0), 22 + oDw (0),11),2}€ + O(e).5 
Here the symbol D is defined by 
2h’ 
D = ————_-- (9.14) 
3(1 — oa?) 


This is a finite quantity; the ordinary flexual rigidity is Dé@E (where E is Young's 
modulus). An example of this type of problem is given below. 

Example. A long rectangular plate is subjected to a uniform tension 7 ()€* on the 
two long edges, and a normal load P?,€? on one face; this normal load does not vary 
along the length of the plate. Find the form of the plate in the strained state. 

In this example, we can neglect the edge effect near the end of the plate by con- 
sidering the plate infinitely long. We assume that the middle surface in the strained 
state is cylindrical, with the generators of the cylinder parallel to the length of the 


plate, that is 
G01 = Woo),11 = Q(x), 022 = Qco)12 = O. (9.15) 


Here x® are rectangular Cartesian coordinates, such that the x?-axis is parallel to the 
long edges, and the x!-axis perpendicular to them; Q is an unknown function. Further- 
more, in this example, 
ae a — so _ yO = 
Po) “ X)fo} 7 Qu) = Qi) X10) = 0. (9.16) 
Then from (9.8b) and the condition that TS are functions of x! only, we have, in 


consequence of the boundary conditions on the two long edges, 


a) = T (2), (9. 17a) ‘2 = Te) =0, (9. 17b) TZ = 0, (9. 17c) 
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Substituting (9.15)—(9.17) into (9.8a), we obtain 


P's) 
Q(x') = ; : (9.18) 
(2) 





Therefore the curvature at any point of the cylindrical surface is proportional to the 
normal pressure at the point. For uniformly distributed pressure, the strained middle 
surface is circular cylindrical. It should be noted that the above conclusion holds in 
general for plates of non-uniform thickness, with the limitation that h is independent 
of x’, 
(b) Type P2: g=0, p=2. FINITE DEFLECTION WITH SMALL EXTENSION IN THE 
MIDDLE SURFACE 


General equations. As in Type P1, we have 
to =n = 3, jo=j =2, ko = k = 2. (9.19) 


By substituting the € series from (8.1), (8.4)—(8.6) into (7.2a, b), (7.4a, b), it is found 
that the lowest power of € occurring in (7.2a, b) is é*, and in (7.4a, b) is e®. The corre- 
sponding coefficients give rise to the equations of equilibrium and compatibility in 
the first approximation as follows: 


70 


— 2ART*@(0)prPi2) ark + BAG" (Qo =a") ton + Pe + 2Xvyoyh 
+ ABT ag (0) 2 (0)28Q (0)0y8® + (Qik) 1x + ie a™qio)nQahk = 0, (9.20a) 
2 Gy (Parah) + Ass "(Qo reQonsh") 1p + $a%*q oprvAlt (Gosh) 1p + P%) 
+ 2Xbyok + (aPa2? + 2a%a*) Qo) 2Qeyrvk + —_ a%*(Qiyh) is = 0, (9.20b) 
mld corarip = 0, (9. 20c); nio] {0} (0) 08 (0) a = 0. (9. 20d) 


The macroscopic tensors in (7.5a, b, c) can be written as 





Ps o ms - 
T?! = {24% Pave + a%Qnh — Ag aweeaonl\e + O(e), (9.21) 
o 


1— 
L? = 3nijarwAtygonshe + O(¢), ° (9.21b) 
T= = {Qiyh + FAT (aeonsh*) ix} 2 + OCC). il 


Here A”, AS” are given as in (6.33a, b). Equations (9.20a, b, c, d) form a set 
of six fundamental equations for the six unknowns pas and qio)as. We see that from 
(9.20d) that the middle surface in the strained state is a developable surface. 

As in Type P1, the problem can be further simplified by introducing wi, such 


that 
J (0)as = Wo) | aB- (9.22) 


We have also 
oC 


a%Qpyh. (9.23) 





TG) = 2AGPayah — AG Qo) 2eQconsh® + "spe 
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We note that 7% is the coefficient of the lowest power of ¢ in T**. It is a symmetric 
tensor, and consequently it has only three independent components. Substituting 
(9.22) and (9.23) into (9.20a, b, c, d), we find that (9.20c) is identically satisfied, while 
the other three equations become 
— TRw yoy + ZAG (W (0) erh*) joy + Pls) + 2X yoy 
a a a 
+ (Qik) |e + a? woo nQak = 0, (9.24a) 
a a 


2 pri r ra 
TE)io + A770) | rwAt) (Woy rah) 1p + Po) + 2X ayo 
a a a a 


+ (a a%? + 2a%a™)w0))2Q(2)7k = 0, (9. 24b) 


Bp a 
1/9] {o]¥ (0) |paW (0) lar = 0. (9. 24c) 
& a 


To sum up, we have for problems of Type P2 a set of four equations (9.24a, b, c), 


in the four unknowns, wy) and 7). 


The special case of uniform thickness will now be treated. Since é is constant, 
(9.24a, b) may be written in the form 


_ Tis (0) 107 + DAAwo) + Pts) + 2X vo 4 + Qinixh + Qf, hAwo) = 0, (9. 25a) 
{TS + Da®ralW (0) |2AW (0) = 2Da*(Awo))?} 19 + Pé + 2XGyo% 

+ (a®Awo) + 2arraP*woyin)Qcarsh = 0, (9.25b) 
and (9.24c) remains unchanged. Here A is the two-dimensional Laplace operator, and 


D is the reduced flexual rigidity as in (9.14). 


Furthermore, when 
(3) = X (2y10) = Qe) = 0, (9.26) 


the equation (9.25b) will be satisfied by putting (¢,) is an arbitrary function of x*) 


Pg Br 1 ~ a 
T%3) = Mio] Mo} P (3) | A —_ Da** a” w (9) |» AW (0) + Da “8 (Aw (o))?. (9.27) 
a a 


And consequently, (9.25a) can be reduced to the form 
— minit}o) 3) | x8 (0) [py + DAAw 0) + $D(Awcoy)* + Pe) 
a a 
+ 2X% igh + O%hdAwo = 0. (9.28) 


Therefore for a plate of uniform thickness under the condition (9.26), we have in 
this type of problem a set of two equations, (9.24c) and (9.28), with two unknowns 
Wo) and gs). In rectangular Cartesians, these two equations may be written as 


W (0),12¥ (0),12 — Wo),11W(0),22 = 0, (9.29a) 
2w 0),120(3),12 — W0),110(3),22 — W0),226(3),11 + DAAwo) 
+ $D(Aw (0)? + Pls) + 2X (21014 + Olahdwo, = 0. (9.29) 


The macroscopic tensors are given by 
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Tis { bcs),22 + $D(wo),22 — W (0,11) AW co) fe? + O(€), 

T? = T2 = — {6,21 + Dw (0) ,12Aw 0) |< + O(4), (9.30a) 
T? = {6 3), + 3D(won — W (0),22) Aw oy fe + O(e), 

LY = —L*® = — D1 — o)we),.2€* + Oe), 

L® = D(wo),11 + ow o),22)€* + O(e), (9.30b) 
L* = — D(wo),22 + ow), 1)e® + O(e), 

T!? = D(Awo)) xe + Oe), T? = D(Aw o)),2€° + OC). (9.30c) 


An interesting example of this type of problem is given below. 

Example: A long rectangular plate of uniform thickness is deformed under the 
actions of (a) uniform tensions 7\3)4€, T3)g€8 and uniform bending moments Lae’, 
Ls)n€ on the two long edges, (b) a normal load P%,)e* on one face (this load does not 
vary along the length of the plate). Assuming that p=2, g=0, find the form of the 
middle surface in the strained state. 

In this example, we can neglect the edge effect at the two ends by considering 
the plate infinitely long. Since the given external force system does not vary along 
the length of the plate, we shall assume that strain and stress are constant along 
this direction. Hence in the first place, the deformed surface is cylindrical, with the 
generators of the cylinder parallel to the length of the plate: 


Jou = Woo), = Q(x'),  Qcoyiz = Woo),12 = Qoy22 = Wio),22 = O. (9.31) 


Here x* are rectangular Cartesians, so that x*-axis is parallel to the long sides and 
x'-axis is perpendicular to them. Q is a function of x', to be determined. 

In the second place, T** is a function of x! only. Since the ends of the plate are 
free from tractions, it follows that T!* and T?* vanish everywhere to the third order: 


Tis) ag Ts) = @. (9.32) 
The component 7" can be written as 
TY = THe + O(e), (9.33) 


where 7\3) is a function of x!, to be determined. 
The problem is to determine two unknowns Q and 743) as functions of x! through 
Eqs. (9.25a, b) under the conditions 


P%) = X(2y0) es X(2)(0) = Qt) 7 Qi) = 0. (9.34) 
Substituting (9.31)—(9.34) into (9.25a, b), we have 
= QT 3) + DQ + P%) = 0, (9.35) (T.3) + DQ’) = 0. (9.36) 


Integration of (9.36) gives 
T(3) + DQ? = C. (9.37) 


Here C is a constant to be determined by the conditions on the long edges. Substitut- 

ing W.o),ae from (9.31) into (9.30b, c), we get 
Ll = —-L® =O(e), L*® = DQe+O0(e), L* = — cDQé + O(e), (9.38) 
T! = DQ, + O(e), T?® = O(e'). (9.39) 
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Then (9.37) becomes 


(Li 2 
ee en 
2D 





(9.40) 


, 


where, by definition, L/j,= DQ. This equation is satisfied everywhere throughout the 


plate. Therefore it is also satisfied at the two long edges, and consequently J 3)4, 
T(3)p, Lysya, Lisyp must satisfy the following relation: 
(Z,3))? (L¢3)p)? 

ee, ae (9.41) 


(3)B _ 


2D 2D 





T (3) 


Therefore we conclude that among 7\3)4, T,s)p, Lys)a, Lys)p only three quantities are 
independent; when any three are given, the fourth can be calculated through (9.41). 


Substituting 73) from (9.37) into (9.35), we obtain 
DQ! — CQ + DQ. = (3) (9.42) 
This is a non-linear differential equation of the second order and third degree in Q. 
When the boundary values of 2 are given (or Lys)4, Lys)z are given), the solution is 


uniquely determined. 
If P3,=0, the problem is identical with the problem of the elastica [8]. For then, 


if we introduce the new variable 6, so that 
Q = 61, (9.43) 


equation (9.42) can be written as 


3D(0,1) — C01 + DO = 0. (9.44) 


The second integral of this equation is 


+D(0,.1)? — C = F cos@. (9.45) 


Equation (9.45) is in the same form as the well known equation for the elastica. The 
constant F can be determined by the boundary conditions on the long edges; 6 is a 
physical quantity which denotes the direction of the tangent to the middle surface 
in the strained state. 

The bending of a rectangular sheet of paper into a cylindrical surface by forces 
and couples applied to two opposite edges may be considered as a problem of the 
above type. There is, however, an edge effect in the neighborhood of the free edges. 


(c) Type P3: g=0, p>2. FINITE DEFLECTION WITH NEGLIGIBLE EXTENSION IN 
THE MIDDLE SURFACE 
General equations. As in type P1, P2, we have 


No = n = 3, jo=j =2, ko =k =2. (9.46) 


By substituting the ¢ series from (8.1), (8.4)—(8.6) into (7.2a, b), (7.4a, b), we find 
that the lowest power in (7.2a, b) is é*, and in (7.4a, b) is e®. The corresponding coeffi- 
cients give rise to the equations of equilibrium and compatibility in the first approxi- 


mation as follows: 





1944] INTRINSIC THEORY OF SHELLS AND PLATES 55 


Atiy ‘(qo xh? Dev + As) rons 0) n8Q (0) 079 (0) rok? + Po) 





. 1 — 2¢ a 
+ (Qlay') = + 2X toyfoyt +5 ” — aq Qiah = 0, (9.47a) 


Alay (@ (0) rw (0) anke® )Ip + §a%"qQ (0) ryAiy “aout, + Ph + 2X Gyoyh 


+ oe a”(Qoyh) ip + (aPa%? + 2a®a™) Qe nQieyk = 0, — (9. 47b) 
—¢ a 


nijToar1a = O, (9.47c) nig) (098A oay = 0. (9.47d) 


The macroscopic tensors in (7.5a, b, c) can be written as 





T# = in a aX%Olh —_ — avert} é + O(e), (9.48a) 
L? = a ae + O(é), (9. 48b) 
Ta? = {3AT“(qonsh*)ix + QGh}e + OC). oy 


Equations (9.46a, b, c, d) form a set of six equations involving only three unknowns 
Q0)as, SO the problem is overdetermined. Let us suppose that qg(o)as can be eliminated 
from these six equations; we get a set of three conditions, which may be written in 
the form 


Vii = Vinh, Qi), P%), X (ayo) = 0, G = 1,2, 3). (9.49) 
Equations (9.49) represent the three necessary conditions on the external force sys- 
tem in order that a plate may undergo finite deflection with negligible extension in 
the middle surface. A special example will be considered as follows. 

Example. Under what circumstances can a portion of a plate of uniform thick- 
ness be bent by normal pressure into a cylindrical surface of finite curvature with 
negligible extension in the middle surface? The normal pressure is assumed to be con- 
stant along the generators of the cylinder. 

In this case, 
X (20) = QO) = Pi = 0. (9.50) 


Let us choose the x?-axis in the direction of the generators of the assumed cylindrical 
surface, and the x!-axis in the perpendicular direction. Then we have as in the ex- 


ample of Types P1, P2 





dio = Q(x’), Qoy12 = Qco22 = 9, (9.51) 
and the equations (9.47a, b) become 
3 - — 20 lias 
Da, + ——— Da’ + P's) +7 — att = 0, (9. 52a) 
2(1 — - 
= D(Q?) 1 -+- aint = 0. (9.52b) 


Integration of (9.52b) gives 


ch 
6 ie aie 9.53 
a + Qi.) ( ) 
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where C; is an integration constant. Substituting 2? from (9.53) into (9.52a), we get 


3 — «)(ch)3!? 
; sa (Ci + Qy)*? 





(chD)*!?{ (Ci + Q%2)) 1/2} ur 


2(1 — «)D'2 
‘ i—2efe\™"... 
+ Ps) + rl 3 D h?!?0 (2) = 0. (9.54) 
— 


This is the required condition to be satisfied by Qi, P%s). 
Let us assume that Q° and P® are of the same order of magnitude; then, since 


P 2) =0, we have 

fo) = 0. (9.55) 
Then the condition (9.54) becomes 
(3 — «)(ch)3/? 
2(1 — o)D'/? 





Ps) = — C3/? = constant. (9.56) 


Furthermore, since the right hand side of (9.53) is constant, the plate is bent into a 
circular cylindrical surface; its curvature is given by 


1 = P 1/3 
haaktes ? Seales. o . (9.57) 
(3 —o)D 





When P?,,=0, we get from (9.57) Q=0. Therefore we conclude that it is impossible 
to bend any portion of a plate of uniform thickness into a cylindrical surface of finite 
curvature with negligible extension in the middle surface, if on that portion of the 
plate the surface force is of the fourth order, and the body force of the third order, 
with respect to the thickness of the plate. 


CONCLUSIONS 


A systematic method of approximation based upon the thinness of the plate has 
been developed in this paper. It is found that thin plate problems may be classified 
into twelve types (P1—P12) according to the relative orders of magnitude of pag, 
as and A. In each case, the problem reduces to the solution of a set of partial differ- 
ential equations, different for different types. These differential equations are given 
in Table I. Furthermore, the principal parts of the macroscopic tensors and the orders 
of magnitude of the external forces for each case are given in Table II. Among these 
twelve types, P1—P3 represent the problems of finite deflection, P4—P8 the problems 
of small deflection, P9—P11 the problems of very small deflection and P12 the prob- 
lems of zero deflection. The problems of finite deflection are discussed in section 9; 
these are new problems, and a simple example for each of these types is solved. The 
problems of small deflections, very small deflection, and zero deflection are familiar; 
the detailed discussion of these types is therefore not necessary. However, we may 
note that the theory of generalized plane stress, the Lagrange-Kirchhoff theory of 
“small” deflection, the von Ka4rm4n theory of “large” deflection and the membrane 
problem can be derived respectively from Types P12, P11, PS, P4. 
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APPENDICES 


(i) TaBLE I.—Table of the equations of equilibrium and compatibility of thin plate problems. 























| | (7.2a) | (7.2b) | (7.4b) | (7.4a) 
Types | q bp — 

| | RARE ECR ae eaey zs | Jaa 
Pi 0 1 ey x : «= x x 
P2 0 2 ££ -4-—-e) gf x x x xX x x 
rPs* | 0 } >2 x x ‘. 2 x eo Fs x x 
P4 | =1 1 x x x x x x 
P5 1 2 x x x x x : 2 x 
P6* =1 2q+1 x x x x x x 
Fy? =1 2q+2 x x eS = £ « x x 
P8* =1 >2qg+2 x x : =. ¢ x x 
P9 2 >1 x x x x x x x x 
P10 >2 2<p<2q 1 x x x x x x 
Pil 22 2q x x x x x .2 x 
P12 | 





8 
IV 

” 

tad 

* 





In this table, the following notation is used: 
The terms occurring in the first equation of equilibrium (7.2a) are 





rn os 2A0y doyP ah, I; = ZAG (Qerh*) joys T3 = Boy ewgrsh® 
- . +o ae i 
I° = Po+ 2Xioh + (O*A) |, I; = 71 we ang” O°A. 


The terms occurring in the second and third equations of equilibrium (7.2b) are 
It = 2AGy*(pah) io, 
a 


Iz = §a%"qzryAly (Qnsh®) |p — AG” (Geugrsh*) |p, 
og 
Iz = P* + 2Xiyh + “ys a(0°A) |p, If = (a®qnal t+ 2a%*qry)Q7h. 
om. a 
The terms occurring in the first equation of compatibility (7.4b) are 


0 8 0 B 
Ji = 2nifo\Mio|Ppr| ass Jz = No} 0} Vor as- 
a 


The term occurring in the second and third equations of compatibility (7.4a) is 
Jai = Ani} Fab |: 
a 


On account of the conditions which hold in the various types of problem, some of 
these terms may be negligible in comparison with others. The table shows by the 
symbol ‘x’ those terms which are to be retained in the first approximation for the 
various types. (The overdetermined problems are denoted by ‘*’.) Thus for example, 
for problems of Type P1, we having the following equations of equilibrium and com- 
patibility in the first approximation: 
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N+i2+ 3 = 0, Ni +i3+ If = 0, J; = 0, Jai = 0. 


These equations are written in terms of the small principal parts instead of in terms 
of the finite coefficients of the lowest power in ¢ (see (9.3a, b, c, d)). 


(ii) TABLE II.—Table of the external force system and the macroscopic tensors for 
various types of thin plate problems. 












































fo [= | 
Types| mo n jo j ko | AS aaeeees ace GE Wa: aR Gaetan 
| | | | ¢ |reerseree] uw [ze] a [rere 
ov fi 2 2 = oa ti« al # fel 2 t** 
| 2 2 1 1 EI 2 a x 3 x 3 z x 
P2 | 3 3 ;i 3 ) -- Jie & wi S35 al 8 fe 2 
P3 | : £4 F.3 3 2 3 | x x| 3 ei 2 oe 2 
| 
PA | q+2| 2 | gti] 1 1 q+1 2| x x} qt+3| x | qt2]x 
q+2 2 q+! | 1 1 q+2 2 |-x x | q+3 x |q@t3 |x x 
Ps | 4/3] 34] 2 2 3] x x| 4 |x| 4 |x x 
P6 | o+3 | 2g+2 | g+2| 2g+1 | 2g+1 | gt+2| 2g+2) x x|qt3| x} gt3]x x 
PT | q+3 | 2g+3 | 2q+2 q+2 gt+2 | 2q+2 | 2qg+3) x x x] qt+3 x |q@+3]}x x 
P8& | q+3 | 2g+3 | 2g+2 | 942 q+2 | 2q+2 | 2g+3) x x | g+3 x |@q@t+t3}x x 
P9 }qt+3|} 3 | gt2] 2 2 | q+2| 3 | x x |q+3|x|@qt3}x x 
P10 | q+3| pti} g+2] p p | qt2| pti] x x|gt3| x|qt3|x x 
Pit | a+3 2g+1 | g+2| 2g | 2q | q+2 | 2g+1| x x | qgt+3 |} x | ¢g+3 |x x 
| 
Pi2| o | p+ o | p | p | @ | pti! x x | © o | 








In this table, the following notation is used: 
The terms occurring in the expression (7.5a) for the membrane stress tensor T° 
are denoted by 





o 
Ti? = 2A3"pah, Ti = — AS ™agugeh, Ti? = —— ah. 
ey 
The term occurring in the expression (7.5b) for the bending moment tensor L** 


is denoted by 
LY) 
LY = 3njawAg qh’. 


The terms occurring in the expression (7.5c) for the shearing stress tensor T*° are 
denoted by 
a na ward 
Ty = Q*hA, Ty = Aq) (Qnsh’) 1+. 
a 


Furthermore, 

mo=order of sum of the normal forces acting on the upper and lower boundary 
surfaces, or order of P®, 

n =order of sum of the tangential forces acting on the upper and lower boundary 
surfaces, or order of P¢, 

jo =order of normal component of body force, or order of X{, 

7 =order of tangential component of body force, or order of Xj, 
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kyo=order of difference of normal forces acting on the upper and lower surfaces, 
or order of Q°, 
k =order of difference of tangential components of forces acting on the upper and 
lower boundary surfaces, or order of Q2, 
t =order of membrane stress tensor T?°, 
u =order of bending moment tensor L*, 
l =order of shearing stress tensor T°. 
This table gives (a) the values of mo, 1, jo, j, Ro, Rk, t, u, 1, (b) the principal terms 
in the expressions for T**, L*8, T*° (denoted by ‘x’). The terms not marked with ‘x’ 
are negligible in comparison with those principal terms. It will be noted that there are 
two lines in the table for P1 and also for P4. This is because, in each case, k may have 
two values. 
For example, in the case of Type P1, we have for T°, L**, 
TH=TPL+TS, Lo= Le, 
while for T@°, 


T°) = T? (if k = 1), 


T®=7T¢°+T72 (if k = 2). 


(To be continued) 
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EFFECT OF A SMALL HOLE ON THE STRESSES IN A 
UNIFORMLY LOADED PLATE* 


BY 
MARTIN GREENSPAN 
National Bureau of Standards 


1. Introduction. The National Bureau of Standards has recently made tests on 
steel columns having perforated cover plates.! Most of the perforations were of so- 
called ovaloid shape, i.e., that of a square with a semi-circle erected on each of two 
opposite sides. The tests on the columns included experimental determinations of the 
distribution of stress in the neighborhood of a perforation, and the results obtained 
aroused interest in the development of a theory for the distribution of stress in a large, 
uniformly loaded plate having a single ovaloid hole. 

In this paper an exact solution to this problem is obtained for a hole having any 
boundary of which the equation can be expressed in the parametric form 


x= pcosB +1 cos 38, y = qsin B — rsin 36. (1) 


The plate is supposed in a state of generalized plane stress, the stress* at points re- 
mote from the hole having the constant nor- 
mal components o,=S;:, ¢,=S,, and the con- 
stant shearing component 72, = Ty. 

Eq. (1) represents a closed curve having 
symmetry about the x-axis and about the 

x y-axis. For certain values of p, g, and r the 
curve is simple, i.e., it does not cross itself. 
By adjustment of the values of p, g, andra 
variety of simple closed curves is obtained, 
including a good approximation to an ovaloid 
The dead Tew seqnnuite Gn wt! and a good approximation to a square with 

ovaloid and the full line the approximate oval. Tounded corners, as well as exact ellipses 

oid of Eqs. (1) and (2). (r=0) of any eccentricity. The approximate 
ovaloid obtained by taking 











Fic. 1. Actual and approximate ovaloids. 


p = 2.063, g=1.108, r= — 0.079, (2) 


is shown compared to the actual ovaloid in Fig. 1. The approximate square obtained 
by taking 
P=q=1, r= — 0.14, (3) 


* Received Nov. 13, 1943. 
1 Ambrose H. Stang and Martin Greenspan, J. Research NBS 28, 669, 687; 29, 279; 30, 15, 177, 411 


(1942-43). 
2 The term stress is used throughout to denote the mean value of the stress over the thickness of the 


plate. 
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is shown in Fig. 2. The sides of the square are y 
parallel to the axes of coordinates. By taking 


p=q=1, 1r=0.14, (4) & 2 


the same square, but with the diagonals parallel 
to the axes of coordinates, is obtained. The 
radius of curvature at the mid-point of the fillet 
is about 0.086 times the length of the side of the 
square. 

2. Curvilinear coordinates. If two sets of 
curves are defined by 


fi(x, y) = a, fo(x, y) = B, (5) Mite aid 


then a pair of values (a, 8) defines the points at ee 
which the corresponding curves (5) intersect, Eqs. (1) and (3). 

and (a, 8) are curvilinear coordinates in the 

x, y-plane. As a special case, the functions of Eq. (5) may be obtained by equating 
real and imaginary parts of both sides of 


w = f(z), (6) 


where w=a+i8 and s=x-++1y. In this case the transformation from the w-plane to 
the z-plane is conformal and the two families of Eq. (5) are orthogonal. The expres- 


sion, 




















dz es ‘ 
~~ = — el, (7) 
dw h 

defines the stretch ratio, 1/h, of the transformation, and gives y, the inclination of the 
curve, 8=constant, to the x-axis. 

In the absence of body forces, the condition that the stresses satisfy the conditions 
of equilibrium is that the normal components, ¢, and gg, and the shearing component, 
Tas, can be derived from a stress function, ¢, by means of the relations® 














, Oo 1 /dd dh? dd OI? 
1 ee 8 da ra, 
— ee a | - 
0a? 2 \0B OB da Ja 
8% 1 (a dh? dg Ah? 
a — * oe va da da a 
and the condition that the expressions (8) satisfy the compatibility conditions is 
ve = (+ )o= (+) n(S+Z)e-0 
dx? ay? da? a6?) \da? ap? 
If a function, F, satisfies Laplace’s equation, 
VF = 0, (10) 


3A. E. H. Love, A treatise on the mathematical theory of elasticity, 4th ed., Cambridge, 1927, p. 91. 
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then F, xF, yF, and p*F=(x?+y?)F satisfy Eq. (9). Functions which satisfy Eq. (10) 
are called harmonic functions, those which satisfy Eq. (9), biharmonic functions. 
3. The coordinate system. The solution of the problem is simplified by the use 
of a coordinate system (a, 8) such that Eq. (1) of the boundary of the hole reduces 
to the form a=ap. Such a system is obtained by writing for Eq. (6) 
z =e” + abe + ac*e™, (11) 
or, separating the real and imaginary parts, 


x 


(e* + abe~*) cos B + ac®e—** cos 38, 


x 
aah s | (12) 
y = (e* — abe~*) sin B — ac*e—** sin 38. | 
For constant a, say ao, Eq. (12) reduces to Eq. (1) for the boundary of the hole, where 
) ] ] y 

p = e% + abe~™, g=e% — abe—%s, r = ac*e—*%, (13) 


From Eqs. (2) and (13) it is easily calculated that for the approximate ovaloid of 
Eqs. (1) and (2), 
et = 1.585, ab = 0.758, ac’ = — 0.314. 


By keeping ab and ac? fixed and varying a and B the appropriate coordinate system 
y § ying » 

















y / 
' ef 
27 / 
2 2.5 
3 T 1 
a” a7 
1.585 
1.25 a 
1.027 
77 1) 
2 3 P 4 
3; 7 
a7 1 47 
3 
57 


Fic. 3. Coordinate system for problem of ovaloid hole. 


for the ovaloid is obtained. This system is shown in Fig. 3. The appropriate systems 
for the approximate square with rounded corners are similarly obtained. Figure 4 
shows the system corresponding to the case of Eq. (3) and Fig. 2, and Fig. 5 shows 
the system corresponding to the case of Eq. (4). 
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Fic. 4, Coordinate system for problem of square hole with rounded corners, sides of 
square parallel to Cartesian axes. 


The coordinates (a, 8) approach polar coordinates (p, @) for large a as follows: 


lim a = log p, lim B = 0. (14) 


The values of h? and its derivatives may be computed as follows. From Eq. (11) 
dz 
vat = e” — abe” — 3ac®e—*”, 
Hence from Eq. (7), 
h-* = e?* + a*h?e—2# + Qarc%e-8* — 2ab cos 28 ) 
+ 6a*bc*e—** cos 28 — 6ac*e—** cos 48, 
and 
1 dh? 


~ — ke — o%*e-** — Zac's * , (15) 
2 da 


— 12a*bc*e-4« cos 28 + 6ac*e—** cos 48), 
Oh? 


op 


Il 





— h*(2ab sin 28 — 6a*bc*e—** sin 28 + 12ac*e—** sin 48). 
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4. The boundary conditions. The statement of the problem may be recapitulated 
as follows. There is given a large plate containing a small hole of the shape given by 
Eq. (1). The edge of the hole is free from stress. The plate is in a state of (generalized) 
plane stress and the components of (mean) stress at points remote from the hole are 


62=S,, 0y=Sy, Tzy=Tzy; or in polar coordinates, 


S,+S, S5.-Sy, ’ : 
= + = ae 26 + Tz, sin 20, 
S.+S, S5:-S 
o% = — sees -— — cos 20 — T., sin 20, (16) 
S,— Sy 


sin 20 + Tz, cos 20. 





> oe ose 
a= 


“ 


The boundary conditions may finally be stated as 





To = Tag = 0, (2 = a), 
S.+S, Ss—S, bs 
oe + —— 28 + Tz, sin 28, 
S.+S, S5:-S (17) 
5 = ——_— — ——— cos 28 — Tay sin 28, 
2 2 
i 
Tap = — — sin 28 + Tz, cos 28, (a = ©), 


The last three of Eq. 17 are obtained by substitution of Eq. (14) into Eq. (16). 
5. The stress function. From the harmonic functions e* sin B and e~“ sin 8 may be 
constructed the biharmonic functions ye* sin 8 and ye~¢ sin 8. From Eq. (12) 


y = e* sin B — abe~* sin B — ac*e~** sin 38. 


Hence 
ye* sin B = 362 — 1¢24 cos 28 — 4ab + 4ab cos 28 + 4ac*e—** cos 48 — fac*e—** cos 28, 
¢ 2 2 2 2 2 2 

— 1 cos 28 — tabe—** + 4abe—** cos 28 + fac*e—4* cos 48 — fac*e—** cos 28. 


nie 


ye-* sin B = 
By dropping the harmonic terms from each of these functions and multiplying by 2 
the two biharmonic functions, 

ob. = e?* + ab cos 28 + ac*e—** cos 48, 

Pb 
are obtained. The biharmonic function 


¢- = ye~* cos B + xe~* sin B = sin 28 — ac*e~** sin 28 


— cos 26 — abe~** — ac*e—** cos 28, 


is obtained in similar fashion. 
The biharmonic function p? may be obtained from Eq. (12): 


p? = x2 + y? = e+ a%%e-2* +4 arc%e-% + 2ab cos 28 
+ 2a*bc®e—4* cos 28 + 2ac*e—** cos 48. 
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Fic. 5. Coordinate system for problem of square hole with rounded corners, 
diagonals parallel to Cartesian axes. 


The non-harmonic stress functions required by this problem are 


1 = 2¢a ed 2abdr ‘as p, 


o2 = — $n, 
and 
os = Pes 
or 
o1 = e?* + a*h2e-** — a2c8e—8« 4 2ab cos 28, 
go. = abe~** + cos 28 + ac*e—** cos 28, 


and 
os = sin 28 — ac*e—4* sin 28. 
In addition, the harmonic stress functions, 


¢7 = 


os = e** cos 28, 4 = € 2% cos 28, $5 =a, e°* sin 28, and ds = e~** sin 28 
will be required. 
The complete stress function may be written 
& = Cidi + Cope + C3hs + Cabs + Cabs + Code + Crbz + Cabs, 


where the C’s are to be adjusted so that the stresses derived from @ meet boundary 


(18) 


conditions (17). Also 
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“ = 2C;(e** — a*b2e-2 + 3a%c%e-8*) — 2C2(abe—** + 2ac*e—** cos 28) 

+ 2C3e?* cos 2B — 2Cye~** cos 28 + Cs + 4Ceac*e—** sin 28 

+ 2C7e?* sin 28 — 2Cse~?* sin 28, (19a) 
rel) 
— = — 4C,ab sin 28 — 2C.(sin 28 + ac*e—** sin 28) 
- — 2C3e?* sin 28 — 2Cye—2* sin 28 + 2C¢(cos 28 — ac*e—** cos 28) 

+ 2C7e?* cos 2B + 2Cse—** cos 28, (19b) 
a% 


—— = 4C,(e?# + a*be~* — 9a*c%e-8*) + 4C2(abe~** + 4ac*e~** cos 2B) 


0a 
+ 4C3e?* cos 2B + 4Cye—?* cos 28 — 16Czac*e—** sin 28 
+ 4C7e?* sin 28 + 4Cge—2 sin 28, (19c) 
0* ; 
36 = — 8C,ab cos 2B — 4C2(cos 28 + ac*e—** cos 28) 
0p? 
— 4C3e2* cos 28 — 4Cye~2* cos 28 — 4C,(sin 28 — ac*e—** sin 28) 
— 4C;e?* sin 28 — 4Cse—** sin 28, (19d) 
0°*o ; ; 
—— = 8C,ac*e—** sin 2B — 4C3e?* sin 28 + 4Cye~*% sin 28 
0adB 


+ 8C,ac*e—** cos 28 + 4C7e?* cos 2B — 4Cse~** cos 28. (19e) 


6. The stresses. Substitution of Eqs. (15) and (19) into Eq. (8) gives the stresses 
in the form 


— = 2C;(Aio + Ai2 cos 28 + Ay, cos 48) 
+ 2C2(Ao + Ao cos 28 + Ao, cos 48 + Age cos 68) 
+ 2C3(As0 + Az2 cos 28 + Ag, cos 48 + A3¢6 cos 68) 
+ 2C4(A4go + Age cos 28 + Agy cos 48 + Aye cos 68) 
+ C5(Aso + Ase cos 2B + As, cos 48) 
— 2C¢6(Age sin 28 + Ag, sin 48 + Age sin 68) 
— 2C7(A7z2 sin 28 + Az sin 48 + Az sin 68) 
— 2Cs(Age sin 28 + Agg sin 48 + Age sin 68), (20) 


— = 2C;(Bio + By cos 28 + By, ccs 48 + Big cos 68) 

+ 2C2(Boo + Bao cos 28 + Ba, cos 48 + Bes cos 68) 

— 2C3(Bso + Baz cos 28 + Bs, cos 48 + Bae cos 68) 

— 2C4( Bao + Baz cos 28 + Bay cos 48 + Bag cos 68) 

— C5(Bso + Bsz cos 28 + Bg, cos 48) 

— 2C.(Bez sin 28 + Beg sin 48 + Bee sin 68) 

+ 2C7(Br sin 28 + Bz, sin 48 + Bz sin 68) 

+ 2Cs(Bsz sin 28 + Bey sin 48 + Bee sin 68), (21) 





1944] 


EFFECT OF SMALL HOLES ON STRESSES IN PLATES 67 


= = 12C;(Dis sin 28 + Dy sin 48 + Dig sin 68) 
— 2C2(Dee2 sin 2B + Dey sin 48 + Dog sin 68) 
+ 2C3(D32 sin 2B + D3, sin 48 + Dge sin 68) 
— 2C,(Da sin 28 + Da, sin 48 + Diag sin 68) 
+ 2C5(Dse sin 28 + Ds, sin 48) 
2C6(Deo + Dez cos 28 + Deg cos 48 + Deg cos 68) 
+ 2C7(Dio + Diz cos 2B + Dz, cos 48 + Dz cos 68) 
2Cs(Dsgo + Dse cos 28 + Dg, cos 48 + Deg cos 68), (22) 


in which 


Ao 


Bao = 


Bs 


By = 
Bu = 


Bso 


= e* + 4a7b? — (24a°c6 + 18a%b2c? — atb*)e—** + 24a*b*c%e—8* — 81atc*e—', 


— 4[abe?= — (3a%bc3 — a*b)e-2* + (Ya%bc® — 3a4b?c3)e** + Yatbc%e—e], 
2(3ac* + a*b? — 6a*b?c®e—4« + Ya%c%e—82), 
2ab — (6a*bc*? — a*b*)e—4* + 30a*bc*%e—*2, 


— 2[e2 — (4ac? — a%*)e-2* + 6(a%c® — a%b%c?)e-8¢ — 18a%c%-102 } 


ab — 8a*bc®e—4** + Ya*dc%e—8e, Aos = — 6a7c8e—*, 

3abe?* — 15a7bc%e-22, Azo = — (e4#* — 15ac* + 3a7b? + 45a%c%e-**), 
abe?* — 9a*bc*e-2, Ase = 3ac?, 

3abe—2* — 3a%bc®e—S, Ag = — [3 — (9ac? — a*b*)e~** — 9a%c%e-8], 
abe-2* + 3a*bc*e—®, Age = — 3ac*e—*2, 

e2* — g*he-2* — 27a7c%e— 8, Ase = — 2a*bc*e—*2, Ass = 6ac*e—?2, 


5 2 [er + (a*b? + 4ac*)e—2* + 6a%c%e—8* + 18a%c%—10], 


— (ab — 2a*bc®e-4* — Ya*bc®e—8*), Age = — 6a7c%e—*, 
e** + 3(a2b? + Sac*) + 45a%c%e-**, Ary = — (abe?* — 9a*bc®e-**), Are = — Sac’, 
3 + (ab? + 9ac*)e—4* — Ya2c%e-8«, Ag, = — (abe—?** + 3a*bc®e—®*), Age = 3ac®e**, 


e** + 4a%h? + (24a%c® + 6a*b*c? + atb*)e—4* — 24a4b*c8e-8* — 8latc!*e—"2, 
— 4[abe** — (3a°bc? — a*b*)e-2* — Ya*bc%e-8* + 18a*bc%e-!*], 

— 2(9ac*? — ab? + 6a*b2c3e-44 — 45a%c%-8), Byg = 12a*bc®e—2¢, 

2ab — (6a*bc*? — a*b*)e—** + 6a*bc®e—*¢, 


2[2(ac? — a%b*)e-2* — 3(4a%c® — a*b2c*)e-8 + Da%c%e-10e], 


< 


ab — 16a*bc*e—4* + 9a*bcSe—*<, Bog = 6(ac®e~** — 2a*c%e—**), 

3abe?* — 15a*bc*e—22, Byp = — (e4* — 15ac? + 3a%b? + 45a%c%e-*2), 
abe?* — 9a*bc3e—22, Bzg = 3ac’, 

3abe—2* — 3a*bc3e—8e, By = — [3 — (9ac* — a*b?)e—** — 9a°c%e-82], 
abe~*« +- 3a*bc*e—*, Bus = — 3ac*e—*, 

e?% — g2h2e-24 — 27a%c%e—-8a, Bsz = — 12a*bc*e—*¢, Bs = 6ac*e*¢, 
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Be = 2[2ac%e 2a + 3(a3b?c3 + 4a7c%)e—S + Qarc%e 10a], 

Bos = — (ab + 4a%bc*e~** — 9a%bc®e-**), Bes = — 6(ac*®e—2* + 2a*c%e-*), 
Bro = e4* + 3(S5ac* + a7?) + 45a%c%e-4*, Brg = — (abe?* — 9a*bc*Fe — 
Se, Bua = — 5 ape + 3a*bc®e—**), Bag = Sac*e 
8a) Dig = — arbc*e-*2, 


3,-9,—10a 
’ 


Deo = e* a“ (2ac* + a°b*)e 2a — 3(2a7c® + a°bh*c)e 6a + 9a*c%e 


4a 
’ 


R .= 3+ (9ac* + a°b?)e 4a _ Qg2c be 


Diz = (10a3bc® + a4b3c?)e—8* — 3atbc%e-!9*, Dig = 2(ac*® + 3a%c%e 


Dey = 4a*bc*e~%, Dog = 3(acte—™* + a’cte—), 
Ds. = e* + 15ac* + 3a7b? + 45a7c%-*, Ds, = — (abe — 9a*bc*e-**), Dsg = — Sac’, 


Daz = 3 + (9ac? + a2b?)e-4* — 9a2cSe-8¢, Day = — (abe~** + 3a*bc®e—%), Dag = 3acke—*, 


Ds2 = ab — 3a*bc*e a Dy = 6ac%e am Deo = 12a*bc*%e wie? 

Dez = — [e — (2ac® + a?b*)e—** — 3(2a7c8 + a*h¢)e—8* — Date re “10a | 

Dea = 4a7°bce ta Deg = — 3(acke—@ = a*c®e—*) , 

Dy. = 3(abe** — 5a*bc*%e—**), Dr = — [et — 3(Sac? — ab?) + 45a°c%e |, 
Dz = abe** — 9a*bc*e > Dre = 3ac’, 

Dso = 3(abe~2* — a*bc®e—*), Dsz = — [3 — (9ac? — a*b*)e—** — 9a°*c%e 8a | 
Dass = abe + 3a*bcre 6a. Dag = — 3ac*e—*, 


Boundary conditions (17) are satisfied by substitution for the C’s in Eqs. (20), 
(21), and (22) of 
f= 10.43). Cle ~KS,=2), C= =- 7s. ) 
— 2(1 — acte-*)C, ab(S2 + Sy) — e740(S, — Sy), 
4(1 — ac®e—*)C, = 4a*bc8e*(S, +S i ~ arm + 3ac*)(S, — Sy), (23) 
[e220 — (ac? — a%b*)e~*0 + (3a7c® + abc?) eS 
— 3a%¢%-10 |(S, + S,) — 2ab(S. — Sy), 


ll 


I 


— 2(1 — acte-*)C;, 


= eT iy, ~— 2(1 + acte~%)Cy = (e%%0 — 3ac*)T ry. 


(1 + acte-40)C, 


The case ac’e~*«° = +1, for which some of the C’s in Eq. (23) are infinite, does not 


correspond to a simple curve for a=a and hence is excluded. 
7. Stresses along the inner boundary. The tangential stress in the boundary 


a= is 
oy = (Splenny 


However, it is simpler to compute it as follows. From Eq. (8), 


Ta + OB 0* 0° 


| 0a? Op? 
Hence from Eq. (19), 


Ca + 8 


73 = 4C,[e + a%b’e-* — 9a%c%e-%« — 2ab cos 28] 
42 
+ 4C,[abe-** — (1 — 3ac%e-**) cos 28] — 4C6[1 + 3acte-*] sin 28. (24) 
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Fic. 6. Ovaloid hole, tension parallel to long axis. Distribution of stress along ovaloid boundary. 
The dashed curve shows the distribution of stress for the case of an elliptical boundary having the 
same ratio of major to minor axis and the same rectified length as the ovaloid boundary. 
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Fic. 7. Ovaloid hole, tension parallel to short axis. Distribution of stress along the ovaloid boundary. 
The dashed curve shows the distribution of stress for the case of an elliptical boundary having the 
sane ratio of major to minor axis and the same rectified length as the ovaloid boundary. 
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Fic. 8. “Square” hole, tension parallel to side. 
The dashed curve shows the distribution of stress for the case of a circular boundary having the same 
rectified length as the “square” boundary. 


For a=ao, 2 =0; hence 


y= (78) a=ay = (Ca -|- oxen 
or from Eq. (24), 


Ot 
— = 4C, [eo + ab%e-20 — 9a%c%e-So — 2ab cos 28] 


5 


1 
> + 4C,|abe 2a — (1 — 3ac%e-4%) cos 28] — 4C,[1 + 3acte—*0 | sin 28, (25) 


in which ho denotes the value of h# for a=ap. 

Substitution into Eq. (25) of ho from Eq. (15), of Ci, C2, and Cs from Eq. (23), 
and replacement of the constants a, b, c, and a, by their values obtained from Eq. (13) 
gives, finally 

[(p? + 6rq) sin? B + (q? + 6rp) cos? B — 6r(p + gq) cos? 28 + 9r? Jo, 

+cte@. 
= (S, + S,)(p? sin? B + gq? cos? B — 9r?) — T.,(p + gq)? Be Babes sin 28 
p+qt2r 

(p? rT q’)(S: T Sy) ne (p + q)*(Sz an Sy) 


rar ; [(p — 3r) sin? B — (q — 3r) cos? BJ. (26) 
q — 2r 
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8. Some special cases. The components of stress at any point in the plate may be 
computed from Eqs. (20), (21), (22), and (23). Of especial interest, however, are the 
values of o,, the tangential stress along the inner boundary, a=ap, at points of which 
the numerically greatest normal and shearing stresses may be expected to occur. 

In this section the values of ¢; are computed and shown for several simple cases. 

Case 1 (Fig. 6). Ovaloid hole, tension parallel to long axis. In this case o; is ob- 
tained from Eq. (26) with S,=7,,=0 and , q, and r as given by Eq. (2). Then 


o1 4.915 — 7.133 cos 28 


S, 3.723 — 2.316 cos 28 + cos 48 





Case 2 (Fig. 7). Ovaloid hole, tension parallel to short axis. Here S,= T,,=0. Then 
ot 1.079 + 7.517 cos 28 
Sy 3.723 — 2.316 cos 28 + cos 48 





Case 3 (Fig. 8). “Square” hole, ten- 
sion parallel to side. Here S,=7.,=0 
and p, g, and r are given by Eq. (3). 
Then 


ot .981 — 2.967 cos 28 





Sa 1.401 + cos 48 


Case 4 (Fig. 9). “Square” hole, ten- 
sion parallel to diagonal. Here S,=T zy 
=0 and p, g, and r are as given by 
Eq. (4). Then 

















ot .981 — 1.606 cos 28 


os 1.401 — cos 48 





In each of Figs. 6, 7, 8, and 9 the 
values of o,/S: or o:/S, are plotted 
along the development of one quadrant 
of the inner boundary of the plate. For 
comparison, there is shown by means of 
the dashed curve in each figure the 
distribution of o:/S: or o:/S, for the 
case of an elliptical boundary having 
the same ratio of major axis to minor 
for the case of a circular boundary having the same axis and the same rectified length as 
rectified length as the “square” boundary. the actual boundary. 





Fic. 9. “Square” hole, tension parallel to diagonal. 
The dashed curve shows the distribution of stress 











SOME PRESENT NONLINEAR PROBLEMS OF THE ELECTRICAL 
AND AERONAUTICAL INDUSTRIES’ 


BY 
ERNEST G. KELLER 


Research Laboratory, Curtiss-Wright Corporation 


1. Introduction. The accelerated growth of research in the field of nonlinearity is 
due to different causes. The general advancement of science requires increasingly 
more precise expressions for the laws of science. Accurate nonlinear equations fre- 
quently depart from the linearized or postulated linear equations which have been 
previously used for approximate results. The quest for perfection and generalization 
and the love of difficult investigations by professional mathematicians play a large 
part in this growth. Another incentive is the increasingly exacting requirements of 
modern manufacturing. These requirements are born of the competitive necessity of 
producing ever improved machines and equipmént in the most economical manner. 
The greatest incentive is necessity. Manufactured equipment and devices must be 
designed to work. 

A nonlinear problem? has been defined as “one which, when formulated mathe- 
matically, reduces to (one or) a system of differential, integral, or integro-differential 
equations such that at least one of the three quantities, a derivative, an integral, or 
a dependent variable, is involved transcendentally or algebraically to a power other 
than the first in at least one equation of the system or in at least one boundary con- 
dition of the system.” Of course, in dealing with applied problems, a physical defini- 
tion independent of all mathematical concepts is preferable, but such is difficult to 
formulate. 

Nonlinear problems resolve themselves into two general types, continuous and dis- 
crete. The first type deals with the behavior of quantities in a field or in at least one 
continuous region of space and, more often than not, reduces mathematically to sys- 
tems of nonlinear partial differential equations. Problems relating primarily to this 
field have been treated by Dr. Theodore von Karman in his Josiah Willard Gibbs 
Memorial lecture.* This paper is both a milestone and a beacon of progress in that 
it is an admirable exposition and inventory of the nonlinear problems of continuous 
fields and at the same time an inspiration and invitation to both the engineer and 
mathematician for further advancement in this difficult field. Among other subjects, 
the von Kaérman lecture treated relaxation oscillations, subharmonic resonance, non- 
linear problems in the theory of elasticity in which the hypotheses of (a) small 
deflections are abandoned, (b) Hooke’s law no longer holds, plasticity, hydrodynam- 

1 A Symposium Address before the American Mathematical Society at Stanford University, April 
24, 1943. Manuscript received Aug. 16, 1943. 

2 E. G. Keller, Analytical methods of solving discrete nonlinear problems in electrical engineering, Trans- 


actions of the American Institute of Electrical Engineers, 60, 1194 (1941). 
3 Theodore von Karman, The engineer grapples with nonlinear problems, Bull. Am. Math. Soc., 46, 


615 (1940). 
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ics; and aerodynamics of (a) ideal fluids, (b) viscous fluids, and (c) compressible 
fluids. The bibliography of the paper contains 178 entries. 

The second type of nonlinear problem is called discrete. Discrete nonlinear prob- 
lems are characterized by the fact that they possess only a finite number of degrees 
of freedom. They are frequently reducible mathematically to systems of nonlinear 
ordinary (total) differential equations or to systems of nonlinear integral equations. 

2. Nature of industrial discrete nonlinear problems. Solutions of nonlinear prob- 
lems in industry are usually in the “small”; i.e., a solution of a system is not required 
for every magnitude whatsoever of the parameters involved. In the solution of such 
problems the greatest single body of theory contributing to nonlinear analysis of 
discrete systems is that which grew out of the attempts of the great French and Ger- 
man mathematicians of the last century to solve the problem of three bodies. While 
their objective, in its complete generality, was never realized, the pure mathematics 
developed (theory of,differential equations, convergence, dominant functions, singu- 
larities, removable singularities, etc.) is today directly applicable in the study of non- 
linear equations of electrical circuits, rotating electrical machines, and various non- 
linear dynamical and aerodynamical devices. The two second largest bodies of theory 
are those of nonlinear integral equations as developed by T. Lalesco‘ and others® and 
the methods of Galerkin® and Ritz along with the modifications of these techniques.’ 

There are at least three salient characteristics of nonlinear engineering problems 
which distinguish them from purely theoretical problems. First, oscillograms, differ- 
ential analyzer solutions,® or other records frequently indicate the nature of the solu- 
tion of the mathematical systems in question. Such mechanical or electrical solutions 
for the same system often differ so much among themselves that there is the risk of 
concluding erroneously that the solution is not unique. (For example, the differential 
equations which yield the two solutions represented in Figs. 3 and 4 also possess 
sinusoidal solutions. Yet the solutions are unique; i.e., the solution in (4) is identical 
with the sinusoidal solution.) Of course, there are systems which do not possess a 
unique solution. In general, even when a solution is unique it may have so many 
manifestations that it is often necessary to integrate the system to determine the 
effect of the various parameters involved. A second characteristic of industrial non- 
linear problems is that frequently the methods of mathematics are not powerful 
enough to yield a complete solution of the problem in sufficiently simple form to be 
usable. Tricks and devices, born of physical concepts, must guide the mathematics 
if a usable solution is to be attained. The mathematics is surely necessary and it is 
just as surely not sufficient. The solution is mathematics plus. A third distinction of 
industrial nonlinear problems is the fact that the derivation of the equations of per- 
formance requires, in addition to a knowledge of mathematics, mathematical physics, 


*V. Volterra, Legons sur les équations intégrales, Gauthier-Villars, Paris, 1913, p. 90. 

5H. Galajikian, Bull. Amer. Math. Soc., 19, 342 (1913); also Ann. of Math., 16, 172 (1915); 
E. Schmidt, Math. Ann., 65, 370 (1908). 

6 A. N. Dinnik, Galerkin's method for determining the critical strengths and frequencies of vibrations, 
Aeronautical Engineering, U.S.S.R., 9, No. 5, 99 (1935). Also W. J. Duncan, Galerkin’s method in mechan- 
ics and differential equations, R&M 1798 (1938). 

7 For additional bibliographies see references 2 and 3 above, also the book, E. G. Keller, Mathematics 
of modern engineering, vol. II, Wiley and Sons, New York, 1942, pp. 303-304. These list a total of 302 
papers; and these papers in turn possess bibliographies. 

8 V. Bush and H. L. Hazen, Integraph of differential equations, J. Franklin Inst., 204, 575 (1927). 
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and engineering, inventive ability in thought. A system of equations may be an in- 
vention of the highest order. It is not always necessary to integrate a system of non- 
linear equations. Often it is necessary only to determine under what conditions the 
physical system is stable. Of course, no single stability criterion exists for nonlinear 
systems such as exists for linear systems. When a solution of a nonlinear system can- 
not be obtained with sufficient rapidity or when it can be obtained but is worthless 
because the time consumed in applying it is too great, it may be possible to obtain the 
information desired by integrating a dominant and a “subordinate” system such that 
the solution of the actual problem is bounded or limited by the solutions of the domi- 
nant and subordinate system. The use of dominant and “subordinate” systems will 
be clear in the following problems. 

3. Some representative discrete nonlinear problems of industry. In this paper a 
number of representative nonlinear systems are treated which illustrate the princi- 
ples enumerated in the last section. These systems are either original, appearing here 
for the first time, or else of very recent date. Some of them pertain to electrical manu- 
facturing, others to aircraft development. Although, as stated above, the derivation 
of the equations of a system is often more important than the solution, none of the 
equations considered are derived here. Some systems are derived in the literature and 
to these references are given. The derivations of the remaining ones can not be given 
for military reasons. These are viewed here merely as hypothetical nonlinear systems. 

1. Nonlinear control circuits. As is well known, the volt-ampere characteristic of 
a nonlinear series circuit (Fig. 1) is represented by the curve in Fig. 2. Such circuits 
have numerous industrial applications due to their rugged mechanical simplicity and 
at the same time their electrical sensitivity. The characteristic in Fig. 2 displays the 
fact that there exists a so-called critical or resonant voltage EH» at which the R.M.S. 
value of the current suddenly increases many fold. For a value of E< Ep (see E sin wt 
in Fig. 1) the current is sinusoidal. For E>» the current has the wave form dis- 


played in Fig. 3. 














Li + 
) ¢ R 2 
= 
<x 
uw 
x +f 
& iE 
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E SIN wd RMS VOLTS 
Fic. 1. Fic. 2. Fic. 3. 


In industrial applications Ep is prescribed. It is required to design a circuit which 
will be sensitive for this prescribed value of EZ». A simple slide rule formula is desired 
which will express Ep as a function of the circuit parameters and of the nonlinear 
reactor employed. The equation of performance for the circuit in Fig. 1 is 


di 1 
Li i) dt + Ri + C idt = — Ecos w(t ~~ to). (1) 
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For the range of interest the current 7 is such that the saturation curve of the reactor 
is single-valued and represented by the equation 


H = ki = x — agx® + asa’. (2) 
With 6=wat, Eqs. (1) and (2) yield 


M “ + R(x — a3x* + asx5) + x. f (x — asx? + asx5)d@ = — Ekcos(@—%), (3) 
C 
where, for a given w, M and x, are constants. The integration of the nonlinear Eq. (3) 
and the development of Z» as a function of the parameters of the physical problem 
are carried out elsewhere and need not be repeated here.® 

2. Nonlinear transmission line phenomena. If a series capacitor is employed in 
the primary side of a transmission line to improve the power factor, curious wave 
forms of voltage and current ensue. The system becomes unstable as far as possessing 
a periodic solution is concerned. This is to be expected since the maximum flux den- 
sity attains a value close to that of the knee of the saturation curve if the transformers 
































— 43.75 =L9x27 - 
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are operating efficiently. The addition of series capacitance is thus likely to create an 
unstable system. In this unstable system, the current and voltage taken on an in- 
definitely large number of wave forms such as shown in Fig. 4. Synchronous motors 
which require sinusoidal applied voltages cannot operate’on currents and voltages of 
the type shown in Fig. 4. 

If the capacitor of the system is shunted by a resistor as indicated in Fig. 5, then 
the equations of performance are 


* E. G. Keller, Resonance theory of non-linear control circuits, J. Franklin Inst., 225, 561 (1938) 
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= + Ry L wi (i ] Yd = FE si 6 -+ 0) 
1 + “fe . a t 4 Sin ’ 
1/1 1 2. I ( 0. 


, : ; ; a (4) 
ioRe = a. f (is — i)d0, (%=41/C), 


= x — a3x° + asa°; H = ai; x = B/Bo; By = dB/dH at H = 0, 


at, = 
or 





zs R (x — a3x* + asx) 


aBy—, + | — + Rill — 3asx* + Sase') 
rd oa 9 


d*x £-Bon dx €-(Ri + R:) 
iste 


2 


= 2 2 Ze 
= [Ri + Zz. [Ez cos (6 +b Ao _ tan! =). 
Re R 


Now R, must have the smallest possible value consistent with stability, since it 
represents a perpetual loss of power. There are ten parameters and two variables. 
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There are infinitely many values of the parameters for which the system is unstable 
and equally as many for which it is stable, i.e., for which the solutions are sinusoidal. 
A convenient slide rule formula is desired giving the above smallest value of R: as 
functions of the other nine parameters of the system. The equations of the system are 
derived and solved elsewhere.!° 

3. Nonlinear springs. It is sufficient to say that, in general, the differential equa- 
tions involving nonlinear springs are integrable by hyperelliptic functions." If damp- 
ing is large a combination of variation of parameters and hyperelliptic functions will 
usually afford sufficient accuracy. 

4. Electric locomotive oscillations. Experience classifies the five oscillatory motions 
of an electric locomotive as pitch, roll, plunge, nose, and rear-end lash. The last two 
are especially important because their pronounced existence in a locomotive produces 
a tendency to derail. Considered superficially, characteristic oscillations of an electric 
locomotive would seem to be similar to those of an ordinary vehicle such as an auto- 
mobile, but experimental data and observation indicate the existence of dangerous 


10 FE. G. Keller, Beat theory of non-linear circuits, J. Franklin Inst., 228, No. 3, 319-337 (1939). 
1 A theory of hyperelliptic functions in usable form is given by F. R. Moulton, On certain expansions 


of elliptic, hyperelliptic and related functions, Am. J. Math., 34, 177-202 (1912). 
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nose and rear-end lash which are not common to an automobile. If the tendency to 
nose exists in an electric locomotive and if the locomotive noses for a speed Vo, then 
it will nose for all speeds greater than Vo. Consequently, nosing is not a resonance 
phenomenon and cannot be avoided by running at a slightly different speed. It might 
be supposed that nosing is due to the coning of the wheels or to the staggering of 
the rails or to a combination of these two possible causes. Such causes, however, 
would produce resonance frequencies for definite discrete values of V instead of in- 
stability for all values of V exceeding Vo. Rails on European railroads are not stag- 
gered and yet electric locomotive nosing still persists. The tendency to nose and the 
violence of the oscillation increase with the weight and power of the locomotive. 

In seeking the source of the phenomenon, consider first an elementary experi- 
ment. Let a miniature set of driving wheels and axle be constructed from two rubber 
paste bottle stoppers and a lead pencil. If the miniature drivers are forced down 
against two rulers as rails, if a torque is applied tending to rotate the wheels, and if 
further in the forward motion slight lateral motion is permissible then an oscillating 
torque will be experienced tending to rotate the axle about a line through the center 
of axle and perpendicular to the plane of the track. The creepage forces between the 
rubber wheels and the rails produce an oscillatory torque. 

The weight of an electric locomotive is so great that it effectively rolls on elastic 
wheels on elastic rails. Making use of this fact and whatever additional facts are nec- 
essary the equations of motion™ can be shown to be 


Mo = 0, 

Nh <+Fi- f= a4(= _ ‘) _~ 24(= - :) -Fi-fi- 24(= - ‘), 
J V V 

Mio + Yi(so — Bin) + dalzo — c& + ban) + dolce — c& + ban) + Béo = 0, 


AE + dac(zo + c& + ben) — roc(zo — c& + ban) + KE 


= — bs(F, + Fo + fi + fe) — 265 Z (Hr +I + He) + 6b5 ff, 





(5) 
Bi — X1bs(zo — Bin) + Aabe(so + cE + ben) + Avbo(so — cE + bon) + Rey = 0, 
: : 2 : 2fds 6fb* . 
CE = — dF; — F.) — ddlfi — fe) — —— (rai — 5x) — —F 
V V 
2frb 
= — (f+ v1 t+ vo) + Fin). 
r 
To the accuracy required, the flange forces are given by 
yi\3 yi\5 yi\"7 
F, = H{—) +7(—) +J fore, 
5; 5; 5; 
(6’) 





y2\? y2\5 ye\" 
neal) lS oe 
6; 61 6; 


12 E, G. Keller, Mathematics of modern engineering, vol. II, p. 72. 
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where the constants H, J, J, j, 1, h are determined from force curves, and 6; and 62 


(6”) 


are lengths shown in Fig. 6. 
The variables ;, ye, ys, ys, and § are eliminated from (5) by means of the relations 


V1 = Yo t+ bs& + dsf, yo = vot dst — def, 
Js = Vo 1 hyé + df, % = Fo + hyé ~~ d,f, 
J = yo t Desk, 


where 65, d3, ds, 4; and hz are lengths defined in reference 12. 

If in (5) Fi= F.=fi=fe=0, then the equations are linear; and the solution can 
be written down at once. This solution is either stable or unstable as indicated by the 
roots of the characteristic equation. The nature of the roots are, of course, a function 
of V, the operating speed of the locomotive. Even if the locomotive is unstable with 
vanishing flange forces, it is stable with non-vanishing flange forces. In this case the 
locomotive is operating roughly and damaging the track needlessly. 

In practical applications, then, it is not necessary to integrate the nonlinear sys- 
tem (5). As a check on the validity of the theory, however, it is necessary to integrate 
the nonlinear system and compare the predicted motion with actual motion as deter- 
mined by runs on a test track. Evidently the solution of (5) for -i= F2=fi=f2=0 
cannot be used as a generating solution for the case of the non-vanishing of the flange 
forces because the stability or unstability of this generating solution is carried over 


into the complete solution. 
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Since nosing and rear end lash are the two motions of most importance, it is suffi- 
cient to resort to elementary means. Consequently, F; and f; [Eqs. (6)] are replaced 
by segments of straight lines as shown in Fig. 6, and the second and sixth equations 
of (5) are solved for yp and ¢ by operational methods (general operational methods 
where both initial charges and initial currents exist must be used.) Since the flange 
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forces are taken as functions with discontinuous slopes, the system of differential 
equations and its characteristic equation change as the flange forces, as functions 
of 1, Yo, Vs, Ys, Change at points a, b, c,d, e, f shown in Fig. 6. The first set of initial 
conditions are chosen by trial and error in such a way that the resulting motion is 
periodic in yo and ¢. The solution for a complete cycle is sufficient. The check of the 
theory is the approximate agreement of computed and test periods. 

5. Nonlinear differential equations of dynamic braking of a synchronous machine. 
The equations of dynamic braking are 


dI (E—TIR) [(rso/s)?+ xaxq] wo I? [x2 + (rso/s)?] 
dt L [(rso/s)? + xaxq] 





s* [(rs0/s)? + waxq][(r50/s)? + xaxg]® 
ds byl? xe + (rso/s)? 
dt ws 


(7) 


where 
2K Prix,(*a — xar)so KPr 
i a) ee a ea) 
JT? JT? 
I being the field current, s the rotor speed, ¢ the time in seconds, all other symbols 
being constant parameters. It is desired to obtain an expression for the time of 
stopping of the rotor as a function of the parameters of the machine. 

The last term in the right member of the first of Eqs. (7) has in all cases a magni- 
tude of approximately ten per cent of its predecessor. Thus a solution as a power series 
in a parameter which vanishes with po is to be expected. Neglecting the term con- 
taining wo in (7) and dividing the first equation by the second, we obtain a solution 
of the resulting equation immediately. However, this solution is an implicit function 
of J and s and such that it is solvable explicitly for either J or s only as a slowly con- 
vergent series. An attempted solution by the method of variation of parameters is 
equally cumbersome. 

It is known from oscillograms, however, that both J and s are monotone decreas- 
ing functions of the time for the interval within which (7) is valid. Change of de- 
pendent variables by 


+ 


[= — +J,e"% and s = sje? 





R 
yields 
dy R 4 R(A? — Ab) [ae + (re*)?|[E/R + Iye-¥|® evt4 
oe a ES Tae KB SE, 
dt L  L{[A2 + (re*)?] ° [42 + (re*)? [42 + (re*)2]8 Ihs¢ (8) 
dz [(E/R) + The”)? [x3 + (re*)?] . 
a, Fain PR AE OL A IIc BR 2 TS ET NR e7* 
ri a [A? + (re*)? }*s2 
where , 
Ao = XarXq, A? = X4Xq, A > Ao. 


The number of revolutions before the rotor of the machine comes to rest is 


1 » 1 e 
N= —{ sdt = — f Soe” *dt. (9) 
2rd 0 2r Jo 
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Now it is sufficient for practical purposes to set an upper limit to N as given by (9) 
provided the upper limit is sufficiently small and provided the results display the 
effect of each parameter of the system. To accomplish this (8) may be replaced by a 
simpler system of equations. Evidently, 


[Az + (re?)?] S (A3 + r)e*, [A2 + (re*)?] < (A? + re? (10) 
for 20. Employing (10) in (7) and integrating, we have 


R(A?- At) ft, 





a i e~**dt 
L L(A? +r%) Jo 
1 (r? a x7) wo ae P > 
$= -f —-+ he | ev*=dt, (11) 
sé, (Az + 7r?)(A?2 +77)? Jo LR 


> 2 * 9 
u(r? + x) ('TE . 
z= — — f FE + Ie" | e~**di, 
52(A? +r*)? 0 R 


where the instantaneous values of y and z as given by (11) are always less than those 
given by the solution of (7) for0<ti< a. 
The system (11) is of the form 


t 
u(t) = (bt) +f K, [?, E, wi(€),---, un(é)|d— (k = 1,2,---,n) 
0 


which is Lalesco’s system of nonlinear integral equations. The solution of this is the 
limit of the sequences 


uy” = ¢x(t), 
t 
us = a) + fi Kalt & x, +> dalle (k= 1,2,--+-,m), 
0 


In the present application ¢:(t) = Rt/L and ¢2(t) =0. For small synchronous machines 
the second approximations u{” and u? give values of y and z such that N in (9) is in 
error by five per cent. The integration in (9) is carried out numerically. Because of 
bearing friction and other decelerating factors not included in (7) the upper limit in 
(9) is finite. 

6. A double-valued nonlinear problem. Consider the integration of the equation 


16 + B6+ w[ki0 + ke tan- k3(6 + a)] = 0. (12) 


This equation was derived ingeniously by W. W. Beman to express an important 
phenomenon in aerodynamics. The quantities J, B, u, ki, ke, ks, and a are all positive 
numbers and the plus or minus sign in (@+a@) is used according as 6<0 or 6>0. 
Evidently, for a particular amplitude of 6, Eq. (12) possesses a periodic solution. 
The period and amplitude of this solution are desired. Eq. (12) in the normal form is 


6 = 6:, 6, == (u/T)[ 0 + ko tan! k3(0 + a) | ~~ B6,/T. (13) 
An integral of (13) for B=0 is 


9 1 
—{—— A @ + a) tan! k3(0 + a) — a. log (1 + &3(0, + ay | 


“3 
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or 
A=+V/c— fH), (14) 


where c=f(@)) and where 4) is the maximum positive displacement for t =f. If (14) is 
used as an equation of change of variable, the method of variation of parameters 
yields 

06; Oc 06; 


m B 
ae [ 10 +e ke tan! k3(0 + a) | =< 61; 
Oc Ot ol I I 


whence 


é= — 2g6;/T, (15) 


28 ¢! 28 ¢° 
c= -- f Oidt +d = -=f 6 d6 + d, 
IJ; I Je, 


0 


or 


where d is an arbitrary constant. From the last equation 


dc 2B 2p a 
— = — —6, = — — (+ Ve — f(0)). 16 
a “the 7 SE vers )) (16) 


To determine the signs in (16) it is evident from (15) that c is a decreasing function 
of the time. Consequently, for 6<0 


dc dc dt dc ae 
: se a Oo 
Thus Eq. (16) is 
nid = + sad /t — f(0) (17) 
dé I 


according as 6;<0 or 6:>0. 
For the integration of (17) it is sufficient to replace 1/c—f(6) by k[c—f(@)] where k 
is determined graphically by 


f ve=FO ao = & fle — sola, 


c=f(@0) or c=f(@s) according as 0.5056) or 0 SO050), and Oo, 0¢, and 6 are 
shown in Fig. 7. The curve in Fig. 7 is the solution (14). With this replacement and 
simple integration 


2Bk 8 
cm eee-H0IT oy — —— f e84/1/(9)d0| (0 $0 S04), 
” (18) 
is 2Bk 0 
C= eH BROT! of — —f e011 (9) d0 (64 < 60"), 
6, 


IA 
S 


where ¢o=f(@0) and c¢ =f(@¢ ). The above values of ¢ are substituted in Eq. (14). The 
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solution is periodic when @ is chosen such that 0j turns out to be equal to 69. The 


period of the motion is then given by 
% 6 
P= 2 f se (19) 
6, 


The numerical integration of (19) presents no difficulty at the limits 0) and @¢ 
since 6, in the vicinity of 0) and @¢ can be replaced by an integrable function fo such 


that the limit of (fo/@:) =1 at 6:.=@0 and (:=0/. 








AS a 



































aa 6, ; oT ORF 
r . : 











Fic. 7 Fic. 8. 


zo de 
7. A nonlinear problem of two oleo-pneumatically coupled masses one of which 
is subject to impact. It can be shown without difficulty that the nonlinear differential 
equations of motion of m; and mz shown in Fig. 8 are 














S$ (S — Sm)3($2 — $1)? 
m3, — W1 — — , <i ey x + esi + f(si) = 0, 
1 (So — oy" 2gc?[A(r) |? 
D 
: (20) 
poS p(S — Sm)*($2 — $1)? 
MS. — NW. + - ~ re=reranie Skee ene ms GO, 
E (Sse — 51) i 2gc?[A (r) |? 
D 
where 
[A(r) ]? = 4} R? — [r; + di(s2e — si) ]?}, f(s1) = (Ko + Ross), (21) 


and »=(1—m), 0<m31. In Eqs. (20) and (21), 51, se, and ¢ are the dependent and 
independent variables, the remaining symbols being constants. 

A solution of (20) is desired for the initial conditions $,(0) = $2(0) =vo. The time ¢ 
is counted from the instant when the lower end of the spring is in contact with a 
fixed horizontal surface. For suitable values of R, r;, and b; a graph of [A (r) ]* is either 


Fig. 9a or 9b. 
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Fic. 9b. 


Even in the case where A (r) =7R the follow- 
ing methods, (a) power series solutions in the 
time, (b) numerical integration, (c) successive 
integrations, (d) Galerkin’s, (e) curvature, and 
(f) expansions in power series in parameters, 
become so laborious that they fail for practical 
purposes. Lord Rayleigh"® has given methods of 
handling differential equations linear in all terms 
except containing a damping coefficient which 
depends on the square of the velocity. The ve- 
locity term is supposed small. In Eq. (20) the 
velocity term is small or large dependent upon 
the stage of the motion. 

Consider the curves shown in Figs. 10, 11, 
and 12. It is evident that the forces relevant to 
these curves can be approximated by the arc 


-B, QoQ; and the small number of secant lines Q;Qs, 
RoR, RiR2, and R2R3. The location of Q;, Ro, 
R,, Re, and R; will become evident from physical 
conditions presently discussed. 
’ For the first interval of motion replace Eq. 
(20) by the equations 
Fic. 10. 
§; + Ko + kos; — al[Ao tf mo(Se aa sy) | -— = B(S2 _ $)*, (22) 
So + y[Ao aa mo(Se —_ s) | = 5(S2 = $1), 
where 
poS p(S — Sm)® poS ; p(S — Sm)? 
=» = a y ae 9 = ’ 
” m 2wyc*r*R? ’ me 2wec*n*R? 


% Lord Rayleigh, Theory of sound, (2nd ed.) vol. I, Macmillan, London, 1894, p. 81. 
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and Ao, mo, Ko, and ko are shown in the figures. The ordinate of R; represents the 
value of (1—£)~-!:2° when the system is at rest under the force of gravity. The points 
Rand R; are located such that no ordinate on the secant lines exceeds the correspond- 
ing ordinate on the arc by more than ten per cent. 

The solution of Eq. (20) is now broken up into two time intervals. We reduce 
Eq. (21) to the normal form by the substitutions 


$y = &: + a, So = & + do, §, = §3, £ = b4, 


where a; and a2 are constants such that no constant term remains in the resulting 
differential equations. Then the equations are 


{i= £3, 
—_— = 
$2 — 64) 
: (23) 
£3 = — (amo + kot: + amot: + B(Es — &s)*, 
£4 = ymoki — ymokte — 5(E4 — &3)*, 
and 
" / 
a, = vig + (aAo — K,)| alyAo = ng |, x ho, 


a, = i moly(g + a@Ay — Ko) — a(yAo — ng)| — ko(yAo — ng) } /moy ko, 


and aA 9=XKp in order that s; may not be positive in its initial motion. That is, the 
origin of time is taken to be the instant at which the upward force of the spring S 
is equal to the downward force due to gas pressure on m. 

The general solution of (23) (with squared terms suppressed) is 


£, = A; sin wt + A» COs wit + Ag SiN wet + A4 COS wol, 

f& = dA; sinwyt+ biAccoswit + beAz sin wet + 02A4 COS wot, 
(24 

£3 = w1A, COS wit — w1A¢ SIN wyl + woA3 COS wet — woA 4 sin Wo, 


4= byw41A 1 COS Wt — byw Ae sin wt + bowoA 3 COS Wel — bowed 4 sin wel, 


wre 


where w; and we are the roots of the characteristic equation and 
f 2 
b, = (amo + ko — wi)/amo, be = (amy + ko — we)/amo. 
The nonlinear terms in (22) are taken into account by the method of variation 
of parameters. Employing (24) as equations of change of variables and remembering 
that (24) satisfies 


0&1 02 
on - §3, wy - £4, 
ot ol 
0&3 P Ok, 
— = — (amo + koi: + amoks, — = ymoki — ymoke 
Ot ot 
we have 
Of Oh . Of. Of . 
— A, + — 4A, +— As +— Ay =0, 
OA 1 OAs OA3 OA 4 = 
(25’) 
dfe Of. . Oke . Oks 
— A, +—4:+— As +— A =0, 


aay dA» ‘Ma Abs: 
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dfs. dés ds. dé; 


hy + Ap + — As + — Ay = Blea — 4, 

big? tk a” ae ai (25/”) 
oz . Ok, . Ok, . Ok, . 

— ht — het ht he - 

OA 1 OA 2 OA3 OA, 


The solutions of (25), after some rather lengthy trigonometric manipulations, are 








5 + d6 : : 6 + 5B ik 
A,=— - — (£4 — £3)? cos wy, Az= - — —— (&, — &3)? sin wy, 
w1(by — B;) wi(b, — di) (26) 
{ Bot (g t,)? t, A on (é £3)? sin wet 
A;=-- &4 — &3)° COS wel, # = — - — €3)~ SIN Wol. 
' wo( be = b;) ' : ‘ wo(be ie by) ‘ . : 


The solution of (26) is obtained with sufficient approximation by using a device com- 
mon in celestial mechanics; i.e., for small values of the time, the A; entering (26) 
through £ and & may be considered constants having the values obtained by the 
solution of (24) for £;= —a,, £= —de, £;=£,=v»9 at t=0. Thus the solution of (26), 
to the accuracy required, is reduced to quadratures. Moreover, since the interval for 
which this solution is valid is small (0 S¢<0.01) the trigonometric functions involved 


may be expanded as power series in ¢ before the quadratures are performed. The 
solution of (25) is 
A; =C; + fib (i= 1,---, 4), (27) 
where {;(0)=0. The substitution of (27) in (24) gives the complete solution for 
0<t<h, where wet; <4 and we>a:. The values of C;=A; as determined above. 

The value of [&(t:) —£3(t1)]? locates the point Q; in Fig. 10. The ordinate of Qe 
is %. The ordinate of R; is given by se(t:)—s:(t;). The ordinate of R, is the value of 
(1—£)-!-2 when the air chamber is decreased to 0.7 of its initial value. 




















“1.20 
12 (i-§) 
iH 
10 $(4) 
9 tT 
8 5 
j 
6 a 
5 
2 - Tan % 
= 3 
: fae 
| Kes 
> A, 
Fic. 11. Fic. 12.* 


* In Fig. 12 the origin should be marked So and the second point of intersection of the straight line 
through the origin and the curve should be marked 5S}. 
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In the solution for the second interval (#;<St<t.) of motion it is sufficient to re- 
place the arcs RiR2R;R, and Q:Q2 by the secant lines R;R, and Q,Q2. The quantity 
($2—$1)? on the interval Q,Q2 may be written 


($e —_ $1)? == Bo + no($e —_ $1), (28) 
and (22) becomes 
(p? + amo + ho)si — amos, = g + aAo — BBo + Bno(S2 — $1) — Ko, 


(29) 
— ymos; + (p? + ymo)s2 = ng — yAo + 5Bo — bmo($2 — $1), 


where any additional constants are shown in the figures. At the new origin of time 
for (29), 51(0) = s2(0) =(0 and $,(0) = $1(t1) =Vj\, $2(0) = $o(t1) = V2. 

The proper determination of the constants 5; and 0}, in the substitution s;=£,+)i, 
$2=&+b2 in (29) yields 


(p? + Bnop + amo + kot: — (Bnop + amo)t, = 0, 
— (dnop + ymo)E1 + (p? + Snoop + ymo)é2 = 0. 


While the characteristic equation of (30) is of the fourth degree, yet its roots are 
widely separated in practical cases and quickly found by Graeffe’s method. 

The values of s:=£1+01, so=£+b2 as given by the solution of (30) do not yield 
the equilibrium positions of m; and mez, because when (s2—$:)? becomes small the rela- 
tion (28) and Eqs. (29) are no longer valid. This is no defect of the solution because 
its purpose is the determination of the maximum accelerations acting on m, and mz». 
These maxima occur in the interval 0SSt. The equilibrium positions of m, and mz 
are determined from static considerations. 

A point of special interest is the determination of the effects of the factor ng upon 
the solution. The above solution is constructed with this in mind. 

The roots of the characteristic equation of (30) have special physical significance. 
In practical cases these are usually one or two pairs of complex roots. If there are 
four complex roots, one pair gives a high frequency oscillation of moderate magnitude 
for m. This is to be avoided. 

If [A(r)]? is given by the graph shown in Fig. 9b the above method is still applica- 
ble. The solution is very sensitive with respect to [A(r)].2 Of course, the intervals 
of solution will exceed two in number, but in each interval the value of (s2—$;)? will 
be given by the ordinates of the arc QoQ: or the secant Q,Q2. 

The most complicated process involved in solving (20) is the solution of a quartic 
equation. 

4. Concluding remarks. The seven problems presented above are representative 
of the nonlinear discrete problems of industry in so far as one nonlinear problem can 
represent a group the members of which differ greatly. No bibliography is given for 
the reason stated in footnote 7. 

Methods of handling industrial nonlinear problems of continuous systems arising 
in industry are reserved for a subsequent paper. 


(30) 
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—NOTES— 


A GEOMETRICAL INTERPRETATION OF THE 
RELAXATION METHOD* 


By J. L. SYNGE (The Ohio State University) 


Let a;;, B; (¢,7=1, 2, - - +, ) be given constants such that a;;=a@j and ) 2) aisxix; 
is a positive definite form. Consider the equations 


> aijx; — By = 0 (¢=1,2,---,m). (1) 
j=1 


The solution is easily expressed as a set of quotients of determinants. However, as n 
increases, the task of calculating the determinants becomes excessively burdensome. 
The relaxation method! provides a set of easy steps by which the solution of (1) is 
approached. The method has been compactly described by Temple.’ 

The purpose of the present note is to give a geometrical description of the relaxa- 
tion method. For the trivial case n= 2 the geometrical description may be displayed 
accurately in a diagram. For n=3 a model may be visualized. For »>3 we pass be- 
yond the region of simple concrete geometrical representation, but in many ways 
geometry in an m-space is closely analogous to geometry in 2-space or 3-space, and 
the geometrical description continues to serve as a general guide to procedure. 

Let us regard x; as rectangular Cartesian coordinates in a Euclidean n-space. Let 
us define 

n n 
H(x) =4 DY ajxix; — DY Bix. (2) 
i,j=l i=l 
The equation H(x)=const. represents a family of ellipsoids E; these ellipsoids have 
a common center, common directions for their principal axes, and common values 
for the ratios of their principal axes. They form, in fact, a family of similar and simi- 
larly situated ellipsoids. 

The equations (1) represent a set of planes (i.e., flats of »—1 dimensions). The 
point of intersection of these planes is the common center G of E. Thus the problem 
of solving (1) is the problem of finding the center of an ellipsoid when its equation 
is given. 

It is important to note that H(x) takes a minimum value at G. H is constant over 
each ellipsoid, and increases steadily as we pass out from G. 

It is not possible to define precisely what procedures are to be regarded as per- 
missible. It is a question of ease of computation. Let us follow Southwell and consider 
an approach to G by steps each of which is parallel to one of the axes of coordinates x;. 

Fig. 1 shows Southwell's procedure. It is a schematic diagram in which the ellip- 


* Received Dec. 29, 1943. 
1R. V. Southwell, Relaxation methods in engineering, Oxford, 1940. 
2 G. Temple, Proc. Roy. Soc. London A, 169, 476-500 (1939). 
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soids are represented by circles. (The ellipsoids can of course be transformed into 
concentric spheres by a linear transformation, which however destroys the orthogo- 
nality of perpendicular lines.) We start with an arbitrary point Po (the zero approxi- 
mation). Let Ey be the ellipsoid which 
passes through Po. Through P» we draw 
a straight line Z parallel to one of the 
coordinate axes. Let Q; be the second 
point in which L cuts Eo. Let P; be the 
middle point of the chord PoQ;. Then P,; 
is the first approximation. 

Since the ellipsoid is a convex surface, 
P, lies inside Ey and so H(P;) <H(P%). 
Moreover it is easy to see that PoQ,; is 
tangent at P, to the ellipsoid Z, which 
passes through P;. Thus, of all points 
on the chord P,Q, the point P; gives 
the smallest value of H. 

The process is repeated, starting 
from P;. The second approximation P, 
is the middle point of a chord P,Q: of Fi, 
drawn parallel to another of the coordi- 
nate axes. In this way we get a sequence 
of points Po, P;, -- - . The success of the 
method depends on the rapidity of the 
convergence of this sequence to G. 

In one important respect the above 
procedure is incompletely defined. When 
we have reached P,,, in which of the directions defined by the coordinate axes are we 
to proceed in order to get P41? There are coordinate axes. Of these one cannot be 
used, viz., that which gave the direction of the step PmsPm. But, of the remaining 


x 





n—1 directions, which should we use? 
Gaskell? has suggested the following plan. Write 


n 


C(x) = >> ayjx; — Bi. (3) 


j=1 


Having reached the point Pn, we calculate the quantities C;(P,). Let Ci(Pm) be the 
greatest of these in absolute value. Then we choose for the step PnPm4: the direction 
of the axis of x;. 

This procedure is called the liquidation of the greatest error, since we obtain 
Ci(Pmsi1) =0. It is interesting to see how this result fits into the geometrical discus- 
sion. The plane C;(x)=0 is the plane through G conjugate to the direction of the 
axis x,. The line PnPm4; is parallel to this axis and tangent at Pm4: to one of the 
ellipsoids, Z,,4:. But the point of contact of a line with an ellipsoid lies on the central 
plane conjugate to the direction of the line. Hence Pmn4: lies on C,(x)=0, i.e., 


Ci(Pm4i) =0. 


3 R. E. Gaskell, Quarterly of Applied Mathematics, 1, 237-249 (1943). 
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But it may well be asked whether the quantities C; themselves possess any deep 
significance. It is true that G satisfies C;(x)=0, but the quantity C;(x) for a general 
point does not represent the perpendicular distance of that point from the plane 
C;(x) =0. This perpendicular distance is 


C, 
p(x) = _ leant _ (4) 


n 1/2 
2 
(x c,) 
j=1 


Should we not liquidate the greatest p; rather than the greatest C;? Or is there a 
better plan than either? 

The following plan is suggested. Having reached the point P,,, we have an option 
on m—1 next points. Each of these points lies on an ellipsoid of the family E. Choose 
that point which lies on the innermost ellipsoid. This is equivalent to saying: Choose 
that point which gives the smallest value to H. 

Now’ for a step in the direction of the axis x; the decrease in H is }C?/a;;. This is 
to be made as great as possible, and so we should pick the direction of the step 
P»Pms1 according to the following rule: Proceed in the direction of the axis of x, where 
C?/ax, is the greatest of the quantities Cj/ai; (t=1, 2, +++, mn). 

Thus C?/a;; is made the criterion rather than Gaskell’s C;. The calculation of the 
former quantities involves slightly more computation, but this may be taken care of 
by making the initial transformation 





xf = (a,;)'/2x,. (5) 
Then 
H =} - Q; jX4X; _ ) 2 Bix., (6) 
where a lk 
a;; = 1, ai; = a;;/(aia;;)""?, B; = B;/(a;)"". (7) 
Now, with 
ee ee (8) 
j=l 


the criterion for the direction of the next displacement is C/? or |C/|, the same 
as Gaskell’s. Moreover the transformation from P, to Pn4; is now extremely simple. 


It is® 
é af me rf. 
P or ey 
° ” ’ EPs =. 
Posst Si,*°*, % —CéE(#), + °°. Be. 


2. E. Gaskell, loc. cit., Eq. (23). 
5 R. E. Gaskell, loc. cit. 
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PROPOSED SYMBOLS FOR THE MODIFIED COSINE 
AND EXPONENTIAL INTEGRALS 


By S. A. SCHELKUNOFF (Bell Telephone Laboratories) 


The standard sine and cosine integrals are defined as follows 


; = sin ¢ < = cos t 
Si z = f —- dt, Cizs= f — dt. 
0 t oe 


The cosine integral has a logarithmic singularity at x =0. Now in problems of electro- 
magnetic radiation x is proportional to the frequency but the impedance functions 
involving Ci x are free from logarithmic singularities at f=0. Thus one expects and 
actually encounters logarithmic functions which cancel the singular parts of the cosine 


integrals. 
For this reason the more suitable function is the following modified cosine integral 


which is an entire function. This function has already been used quite frequently, and 
we wish only to suggest that a standard notation be adopted for it. 

Inasmuch as one is frequently interested in the analytic properties of impedance 
functions over the entire oscillation constant plane, the following modified exponen- 


mn z 1 — ev 
Ein z = — — dw. 
0 w 


The independent variable z will be proportional to p=£+iw where w= 2z times the 
frequency. Then, on the imaginary axis we have 


Ein (iy) = Cin y + 7 Si y, 


tial integral is suggested 





where y is proportional to the frequency. 
The even part of Ein z may be designated as Cinh z and the odd part Sih z. 
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BOOK REVIEWS 


The mathematics of physics and chemistry. By Henry Margenau and George Moseley 
Murphy. D. Van Nostrand Company, Inc. New York, 1943. xii+581 pp. $6.50. 


Contents: 1, The mathematics of thermodynamics. 2. Ordinary differential equations. 3. Special 
functions. 4. Vector analysis. 5. Vectors and curvilinear coordinates. 6. Calculus of variations. 7. Partial 
differential equations of classical physics. 8. Eigenvalues and eigenfunctions. 9. Mechanics of molecules. 
10. Matrices and matrix algebra. 11. Quantum mechanics. 12. Statistical mechanics. 13. Numerical cal- 
culations. 14. Linear integral equations. 15. Group theory. 

The need for comprehensive manuals of mathematical tools is widely felt by workers in various ap- 
plied fields. The readers of this Quarterly may, therefore, envy the theoretical physicists and chemists 
for whom the present book is primarily intended. However, it appears from the above table of contents 
that the book covers such a great variety of topics that almost everyone will find some chapter of par- 
ticular interest. In this connection the chapters on special functions and special coordinate systems de- 
serve particular mention. 

The authors have well succeeeded in making the book appear as a homogeneous unit although the 
individual chapters are independent and show a refreshing lack of formal uniformity. In some chapters 
physical theories are treated at very considerable length, while other chapters are quite mathematical in 
form. Formal deductions are given in general, but often it seemed more desirable merely to record formulas 
or facts. “The degree of difficulty of the treatment is such that a Senior majoring in physics or chemistry 
would be able to read most parts of the book with understanding.” 

Occasionally, a more daring departure from customary lines would have made the book still more 
useful. Thus some numerical methods which are often presented and hardly ever used would better have 
been omitted in favor of a more thorough presentation of the really useful techniques. The modern statis- 
tician will regret to find the theory of errors treated along conventional, obsolete lines. The magic spell 
of purely conventional but impressive terms such as “probable error” has proved very dangerous indeed 
and inspires an unjustified confidence. The physicist who still believes in the normalcy of observational 
errors should consult W. A. Shewhart’s “Statistical Method From the Viewpoint of Quality Control” 
(Washington 1939). There, starting on p. 66, he will find a most interesting analysis of some measurements 
among the very elite (velocity of light, the gravitational constant, Planck's constant). They all show com- 
plete lack of statistical control, and even the simplest methods of industrial quality control could be used 
for an improvement. 

In general, the presentation is very clear. Only occasionally an attempt at mathematical scphistica- 
tion makes itself felt. Thus the authors first introduce vectors in the usual (most satisfactory) manner. 
Then (pp. 134-135), rather unclear references are made to a more restrictive analytical definition. The 
passage culminates in the puzzling statement that [y, x] (which, by the way, is the gradient of the func- 
tion xy) “does not define a vector.” It does. And the authors themselves make free use of gradients and, 
on the other hand, they (p. 135) “do assume that all of the vectors discussed are proper vectors.” 

W. FELLER 


Navigational trigonometry. By P. R. Rider and Ch. A. Hutchinson. The Macmillan 
Company. New York, 1943. ix+232 pp. $2.00. 


The reviewer has considered this book more from the standpoint of a person studying the principles 
underlying the art of navigation, either for the first time or as a refresher, than as a mathematical text- 
book. 

The book, as the authors say, is “a revision and expansion of part of Rider’s Plane and Spherical 
Trigonometry.” The general arrangement of the material is very good, both as to the sequence of topics 
taken up by chapters and the presentation of the material in each chapter itself. Chapter by chapter it 
leads the student from fundamental definitions through the solutions of right spherical triangles and 
oblique spherical triangles which are necessary for the student to know if he is to thoroughly understand 
his navigation. Admitting that one can learn to navigate and use the short cuts common to practical 
navigation without a very thorough background of spherical trigonometry, nevertheless the more com- 
plete his knowledge of this branch of mathematics, the better navigator he will be and the more he will 
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enjoy working out navigational problems. This phase has, in the reviewer's opinion, been very well 
handled by the authors, who have shown good judgment in maintaining the proper balance between the 
amount of detail used in “proofs” and the confidence shown in the intelligence of the student in assuming 
that he will either accept certain facts or will be able to complete the detailed proofs himself. 

The chapters on The Terrestrial Sphere, Charts, The Sailings, Astronomical Triangle and Lines of 
Position are presented in clear, concise English and in logical order, giving the student the information 
necessary for him to understand the problems which will confront him later when he takes up navigation 
as a working tool. The authors very sensibly do not attempt to include in these chapters everything that 
a man must know in order to actually navigate, but leave that to other books written especially for this 
purpose. 

Throughout the book the method of presentation of material is excellent. Each chapter contains cer- 
tain proofs and facts followed by problems or exercises based on preceding information, giving the student 
an opportunity to apply the principles discussed. The fact that the answers to certain problems are given 
in the back of the book gives the student the chance to know whether or not he has used the proper 
method of solution, and also the satisfaction of knowing that he has successfully accomplished his task. 

The inclusion of an appendix discussing briefly the standards of accuracy is, in the reviewer's opinion, 
very well worth while. This subject, often neglected, is not well understood by students who have had 
little experience in mathematics, and is all too often not recognized even by those who have had such 
experience. 

The problems throughout the book are well thought out and the authors have given careful study 
to the matter, laying special emphasis on the authenticity of materials and assumptions so that the prob- 
lems are as practical as possible in a text of this size. 

The book contains a complete five place table of natural and logarithmic haversines with one minute 
intervals, which is a notable and welcome innovation in a textbook on trigonometry. It has been the re- 
viewer’s experience that the beginner finds it confusing to use the table in Bowditch with its variable 
interval. This table, along with the table of Common Logarithms of the Trigonometric Functions makes 
it necessary for the student to make less frequent use of Bowditch, which, because of its size, is rather 
awkward to manipulate. 

All in all, the authors have accomplished what they set out to do. The book fulfills their claims even 
better than might be expected and should prove to be very popular in the teaching and studying of the 


basic mathematical problems underlying the principles of navigation. 
LEIGHTON T. BoHL 


Table of circular and hyperbolic tangents and cotangents for radian arguments. Prepared 
by the Mathematical Tables Project, Work Projects Administration of the Fed- 
eral Works Agency; conducted under the sponsorship of the National Bureau of 
Standards. Official Sponsor: Lyman J. Briggs. Technical Director: Arnold N. 
Lowan. Columbia University Press. New York, 1943. xxxviii+410 pp. $5.00. 


The main table gives the values of tan x, tanh x, cot x and coth x over the range x=0 to x=2 at in- 
tervals of 0.0001. Circular and hyperbolic tangents are given to 8 significant figures for 0<x $0.01 and 
for 0.1 <x <2, and to 9 decimal places for 0.01 <x <0.1. Circular cotangents are given to 8 significant 
figures for 0.1 $x 1.57 and 1.575 Sx $2, to 8 decimal places for 0<x 0.1 and to 13 decimal places for 
1.57 <x <1.575. Hyperbolic cotangents are given to 8 decimal places for 0 $x $0.1 and to 8 significant 
figures for 0.1 Sx <2. The second central differences for all these functions are given wherever linear inter- 
polation is not sufficient. Auxiliary tables contain the values of the circular and hyperbolic tangents and 
cotangents to 10 decimal places over the range x =0 to x=10 at intervals of 0.1; the values of the inter- 
polation coefficients for the formulas of Gregory-Newton and of Everett; the values of m/2 for integer 
values of m from 1 to 100; and values facilitating the conversion from radians to degrees and from de- 


grees to radians. 
W. PRAGER 


























